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The  views  and  conclusions  contained  in  this  document  are  those 
of  the  authors  and  should  not  be  interpreted  as  necessarily  representing 
tie  official  policies,  either  expressed  or  implied,  of  the  Advanced 
Research  Projects  Agency  of  the  U.  S.  Government. 


ABSTRACT 


A model  for  the  prediction  of  the  temperature  and  vapor  fields 
created  about  a small  water  droplet  undergoing  irradiation  by  a laser 
beam  has  been  developed.  Time  dependent  and  steady  state  solutions 
ot  the  model  are  discussed,  as  well  as  the  possible  effects  of  a number 
of  aerosol  properties  not  included  in  the  model.  Estimates  of  character- 
istic phase  shifts  to  be  expected  in  propagating  through  standard  atmospheric 
aerosol  distributions  are  also  presented.  While  the  model  is  quite  general, 
the  calculations  in  this  report  are  limited  to  DF  laser  wavelengths. 
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I.  INTRODUCTION 


There  is  considerable  current  interest  in  the  effect  of 
atmospheric  aerosols  on  laser  propagation1"4.  In  the  case  of 
propagation  of  low  intensity  radiation  through  the  atmosphere  the 
primary  extinction  mechanisms  are  absorption  and  scattering  by 
both  aerosols  and  molecules.  As  the  radiation  intensity  increases 
the  interaction  between  these  extinction  mechanisms  and  the  radiation 
field  can  become  important.  The  phenomenon  of  thermal  blooming 
occurs  when  the  laser  energy  absorbed  by  atmospheric  aerosols  or 
molecules  is  converted  into  translational  energy  of  the  surrounding 
air.  The  index  of  refraction  of  this  heated  air  will  differ  from  ambient 
which  in  turn  can  cause  beam  de-focusing  and/or  steering.  The  magnitude 
of  this  effect  will  depend  upon  both  the  laser  pulse  shape  and  the  heating 
phenomenology. 

The  case  of  laser  blooming  resulting  from  molecular  absorption 
has  been  treated  in  some  detail  and  is  included  in  the  larger  laser 
propagation  codes  (for  example,  Ref.  6).  The  atmospheric  heating  due 
to  molecular  absorption  may  be  treated  relatively  straightforwardly  in  that 
it  occurs  continuously  and  uniformly  in  space  and  the  temperature  increase 
at  any  one  point  is  directly  proportional  to  the  local  laser  intensity.  Recently, 
the  heating  resulting  from  the  absorption  of  laser  energy  by  aerosols  has 
received  some  attention.  U3’  7 The  situation  in  this  case  is  more  com- 
plicated than  that  of  molecular  absorption  in  that  the  absorption  and 
resulting  heating  occurs  at  discrete  points  rather  than  continuously  and  thus 
the  temperature  fields  in  the  surrounding  air  are  highly  nonuniform. 
Blooming/scattering  predictions  have  been  presented1"3  for  the  temperature 
fields  resulting  from  the  absorption  of  laser  radiation  by  nonvaporizing 
aerosols.  As  will  be  discussed  further  in  the  text  the  main  features  of 
these  heating  models  are  that  the  rate  of  temperature  increase  of  a given 
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aerosol  is  linearly  proportional  to  the  laser  intensity  and  that  the 
increase  in  temperature  above  ambient  about  the  aerosol  will  scale 
inversely  with  distance  from  the  aerosol  out  to  some  distance  which  is 
an  increasing  function  of  time. 

The  temperature  fields  resulting  from  the  absorption  of  laser 
energy  by  a vaporizable  aerosol  have  been  discussed  briefly  in  Ref.  7 
and  will  be  considered  in  more  detail  here.  As  discussed  earlier'  the 
temperature  increase  of  the  aerosol  need  no  longer  scale  linearly  with 
laser  intensity  and,  furthermore,  the  temperature  increase  in  the  air 
around  the  aerosol  will  not  necessarily  decrease  invers  ely  with  radius. 

The  subsequent  prediction  of  blooming/scattering  resulting  from  such 
temperature  fields  necessarily  becomes  more  complex  than  for  the  case 
of  nonvaporizing  aerosols. 

This  report  is  concerned  with  the  prediction  of  the  temperature 
and  vapor  fields  arising  from  the  interaction  of  a laser  beam  with  hygro- 
scopic aerosols  at  sea  level.  Although  the  modeling  presented  is  quite 
general  the  discussion  will  be  limited  to  DF  laser  wavelengths.  Section  II 
provides  a brief  description  of  atmospheric  aerosols  and  their  optical 
properties.  The  details  of  the  mathematical  model  describing  the  inter- 
action of  a laser  beam  w'ith  an  aerosol  (taken  as  pure  I-^O)  are  presented 
in  Section  III.  Simplified  analytic  solutions  to  the  model  as  well  as 
detailed  computer  predictions  for  a matrix  of  laser  intensities  and 
aerosol  sizes  are  given  in  Section  IV.  A number  of  real  aerosol  effects 
not  included  in  the  modeling  are  discussed  in  Section  V.  These  include 
nonuniform  aerosol  heating,  original  solute  effects  and  ebullition. 

Section  VI  provides  a brief  discussion  on  the  effect  of  the  laser-aerosol 
interaction  on  laser  propagation.  Lastly,  a summary  of  the  report  is 
presented  in  Section  VII. 


II.  PROPERTIES  OF  ATMOSPHERIC  AEROSOLS 


A-  Aerosol  Particle  Size  Distributions 

The  particle  size  distribution,  mass  loading  and  composition  of 
natural  sea  level  aerosols  vary  considerably  depending  upon  such  factors 
as  time  of  day.  geographical  location,  weather,  etc.8  The  standard 
continental  distribution  ’ will  be  used  in  this  section  for  computational 
purposes  because  of  the  simplicity  of  its  mathematical  form.  This 
distribution  is  defined  by  the  relationship 
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where  N'A  is  the  number  of  particles/cc  of  radius  R.  The  R*4  dependence 
of  the  distribution  has  been  verified  by  a number  of  measurements’0: 
however,  the  functional  form  can  be  somewhat  misleading.  Assuming 
an  average  aerosol  material  density  of  1.5  gms/ce  the  total  mass  of  the 
aerosol  distribution  represented  by  Eq.  (1)  iS  a modest  12  ~gm/m3. 
furthermore,  the  distribution  corresponds  to  a total  particle  density  of 
1 10  part/cc  110'  cm  < R < •,  x 10'3  cm.  ) 

Oceanic  aerosol  distributions  are  not  as  well  defined  as  the 
continental;  however,  there  is  a 'standard'  maritime  distri  tuition4,  “ 
which,  although  somewhat  different  in  shape  than  the  continental  distribution, 
also  cirreep  nds  to  a mass  loading  of  approximately  12  -gm/m3.  There 
can  be  a large  number  of  particles  with  radius  < to'°  cm  in  both  the 
continental  and  maritime  environments;  however,  for  most  situations 
these  will  not  constitute  a significant  fraction  of  the  total  aerosol  mass 
loading  and  have  been  neglected  in  the  present  analysis. 
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The  variation  in  aerosol  mass  loading  from  the  ' standard" 
distributions  can  be  quite  large.  As  examples,  aerosol  mass  loadings 
in  urban  areas  such  as  Los  Angeles,  are  typically  measured  in  hundreds 
of-gm  m and  sea  salt  concentrations  in  Pacific  Ocean  aerosols  have 
been  found  to  sary  from  10  - 1000  h gm/m3  depending  on  wind  velocity.  ‘ 1 
Of  course,  scenarios  involving  mist,  fog  or  even  ship  induced  spray  could 
result  in  still  higher  mass  loadings. 


B.  Relative  Humidity  Effects 

it  has  been  demonstrated  that  hygroscopic  aerosols  increase  in 
size  with  increasing  relative  humidity,  this  accretion  resulting  from  the 
condensation  uf  water  vapor  on  the  aerosol.  A relationship  widely  used4,  9 
to  predict  the  factor  F by  which  the  radius  of  the  aerosol  increases  with 
increasing  relative  humidity  is 


F = 1 - 0.  o In  (1  - R.H.  / 100) 


(2) 


’.'.here  R.H.  is  relative  humidity.  This  relationship  is  obviously  unsatis- 
factory for  R.H.  = 100;  however,  it  is  expected  to  provide  reasonable 
estimates  of  aerosol  growth  at  high  relative  humidities.  Unfortunately 

there  is  little  experimental  data  available  to  provide  validation  for  Eq.  (2) 

12 

Recently  DeBary  et  al  presented  measurements  on  the  growth  of  sea 
spray  aerosols  with  relative  humidity.  No  details  of  the  measurement 
technique  were  given.  These  measurements  are  shown  in  Fig.  1 along  with 
the  growth  curve  predicted  by  Eq.  (2).  As  can  be  seen  the  data  fall  some- 
what below  the  prediction  over  the  full  range  of  relative  humidity.  This 
difference  corresponds  to  approximately  a factor  of  two  in  aerosol  volume 
at  *8  relative  humidity.  Note  there  is  a hysteresis  effect  in  the  aerosol 
growth  measurements  at  relative  humidities  less  than  80%  (i.e.  the  growth 
parameter  varies  depending  upon  whether  the  measurement  was  made  with 
increasing  or  decreasing  relative  humidity). 


I0'«  10° 


(I-R.H./I00) 

Growth  Factor  for  Hygroscopic  Aerosols  vs.  Relative  Humidity 
, Data  for  sea  spray  from  Ref.  12,  - Eq.  (2)  of  text. 


It  should  be  noted  that  there  is  some  theoretical  justification  for 
the  use  of  Eq.  (2).  A pure  water  droplet  cannot  attain  an  equilibrium  con- 
figuration unless  the  surrounding  water  vapor  concentration  corresponds  to 
super-saturation.  On  the  other  hand  salt  solution  aerosols  can  reach  equili- 
brium configurations  at  relative  humidities  below  100%  because  of  the 
reduction  in  the  water  vapor  pressure  at  the  droplet  surface  caused  by  the 
presence  of  the  salt.  A salt  solution  aerosol  will  reach  a saturation  level 
at  a relative  humidity  somewhere  below  100%.  As  the  relative  humidity  is 
increased  above  this  level  the  solution  will  be  diluted  and  the  particle  will 
grow  to  a new  equilibrium  size.  The  resulting  growth  law  will  be  a function 
of  the  solution  properties  and  for  the  case  of  NaCl  this  growth  is  in  reasonable 
agreement  with  the  predictions  of  Eq.  (2)  at  relative  humidities  above  80%, 
this  latter  condition  corresponding  to  the  saturated  solution.  The  details 
of  this  phenomenon  are  given  in  Section  V. 

In  any  case  this  relationship  only  pertains  to  hygroscopic  aerosols, 
and  the  percentage  of  aerosol  mass  loading  corresponding  to  hygroscopic 
materials  can  vary  significantly.  Maritime  aerosols  can  be  dominantly 
hygroscopic  whereas  aerosols  over  urban  areas  may  be  composed  largely 
of  non-hygroscopic  materials.  (For  example  it  has  been  reported  that  only 
15  o of  the  aerosol  mass  loading  over  Los  Angeles  is  hygroscopic13).  The 
remainder  of  this  discussion  will  be  limited  to  hygroscopic  aerosols.  It 
is  clear  for  this  case  that  relationship  (2)  provides  for  a dramatic  variation 
in  aerosol  mass  density  (or  volume)  with  relative  humidity.  For  example 
at  98%  relative  humidity  the  total  aerosol  mass  density  for  the  continental 
distribution  would  be  740  ilgm/m3  as  compared  to  12  Hgm/m3  on  a dry 
day.  It  should  be  emphasized  that  this  apparently  high  mass  loading  still 
corresponds  to  a relatively  clear  day,  the  scattering  coefficient  for  visible 
light  being  < 10  1 km"1  . 

C.  Particle  Absorption  Cross  Sections 

If  the  aerosol  particles  are  assumed  to  be  spheres  their  absorption 
cross  section-  n ay  be  calculated  by  use  of  Mie  Theory  (particle  radius  of 
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same  order  as  laser  wavelength.  ).  Thus  the  absorption  cross  section  is 
a complicated  function  of  the  complex  index  of  refraction  of  the  absorbing 
material,  m = n - in',  and  the  ratio  of  particle  size  to  wavelength,  R/x  . 

The  simplest  case  corresponds  to  when  | m - 1 J <<  1 and 

4TTn'  R/X  < 1 whence  it  can  be  shown14  that  the  absorption  cross  section 
is  given  by 

Q = 4/3  IT  r3  (3) 

with  the  quantity  4tt  n’/X  defined  as  & , the  bulk  absorption  coefficient. 

This  limit  is  generally  referred  to  as  volume  absorption.  Unfortunately, 
real  aerosol  materials  exhibit  values  of  n ranging  from  1.  34  to  1.  5 in  the 
infrared  and  prediction  of  the  absorption  cross  section  requires  use  of 
the  full  Mie  theory. 

Suc.i  predictions  have  been  made  for  the  case  of  water  droplets, 

with  in  - 1.  364  - 0.  0034i,  for  a wavelength  of  3.  8 M-m  (corresponding  to 

a strongly  transmitting  DF  laser  line).  These  predictions,  normalized 

by  particle  volume,  are  shown  vs.  aerosol  radius  in  Fig.  2.  Shown  for 

comparison  are  the  normalized  absorption  cross  sections  resulting  from 

volume  absorption,  Eq.  (3),  and  that  predicted  by  geometric  optics,  which 

of  course  is  only  valid  when  X <<  R.  As  can  be  seen,  over  much  of  the 

aerosol  size  range  of  interest  the  Mie  theory  prediction  is  larger  than  that 

due  to  volume  absorption.  At  larger  particle  sizeT  the  absorption  coefficient 

falls  off  because  of  self  absorption  and  at  the  smallest  radii  asymptotes  to 

14 

the  Rayleigh  limit  , i.  e. 


Q ~ - 


(4) 
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At  approxi  'ately  R = 60  pm  the  Mie  theory  and  geometric  optics 
predictions  converge.  It  should  be  noted  that  the  geometric  prediction 
rises  above  that  of  volume  absorption  because  of  nonuniforrn  absorption 
within  the  aerosol  resulting  from  focusing  effects.  It  may  be  that  the 
increase  in  the  Mie  theory  absorption  cross  section  over  that  of  volume 
heating  may  also  result  from  nonuniform  absorption  within  the  aerosol. 

(This  can  only  be  verified  by  an  examination  of  the  electro  magnetic  field 
within  the  aerosol  as  predicted  by  Mie  theory).  In  any  event  in  the  modeling 
to  be  presented  in  Sections  III  and  IV  uniform  absorption  has  been  assumed 
for  simplicity.  A discussion  of  possible  effects  arising  from  nonuniform 
absorption  is  presented  in  Section  V. 


D.  Atmospheric  Aerosol  Absorption  Coefficients 


The  total  atmospheric  aerosol  absorption  coefficient  is  defined  by 
R„ 


/ 


k - / Q (R,  n , n)  (dN  /dR)  dR 


(5) 


R 


1 


where  dN^/dR  is  now  the  aerosol  distribution  function,  Eq.  (1),  with  the 

effects  of  relative  humidity  included,  (Eq.  (2)  ).  For  simplicity  the  aerosol 

has  been  taken  to  be  composed  completely  of  hygroscopic  material  havin 

the  complex  index  of  refraction  of  water.  This  approximation  is  most  relevant 

to  a marine  environment  on  a humid  day.  (It  has  been  previously  denion- 
4 

strated  that  the  "standard"  continental  and  maritime  distributions  result 

in  similar  absorption  coefficients  for  a.  4 pm). 

There  is  still  some  uncertainty  in  the  imaginary  part  of  the  index  of 

refraction  of  liquid  water  at  3.  8 pm.  In  two  recent  reviews  the  values  of 
-3  . - 3 

3.4  x 10  and  6 x 10  have  been  recommended  (Refs.  1,5  and  16,  respectively). 
Absorption  coefficients  for  both  values  of  n'  have  been  calculated  in  order 
to  bound  the  actual  aerosol  absorption  coefficient.  For  simplicity  the 
calculations  were  performed  assuming  volume  heating,  Eq.  (3),  and  the 
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prediction  tor  = 6 x ID’3  was  increased  hy  a factor  of  i.  5 as  a crude 

correction  factor  for  Eq.  (3)  to  account  for  the  fact  that  n,  - 1.36. 

These  predictions  are  shown  vs.  reli  five  humidity  in  Fig.  3 along 
with  similar  predictions  for  an  aerosol  distribution  with  ten  times  the 
mass  loading.  The  shaded  area  between  each  set  of  curves  approximately 
represents  the  uncertainty  in  the  absorption  coefficient.  It  is  to  be  noted 
that  these  predictions  are  much  lower  than  those  presented  by  Hodges  4 

for  the  same  aerosol  distribution  function.  The  cause  for  this  discrepancy 
is  not  clear. 

An  important  consideration  is  the  determination  of  the  conditions 
under  "hich  absorption  will  dominate  molecular  absorption.  McClatchey 
and  Selby  have  made  predictions  of  atmospheric  laser  attenuation  for  a 
number  of  important  HF  and  DF  laser  lines.  An  example  of  a strong  DF 
User  llne  with  weak  molecular  absorption  is  the  v = 2-  1,  J = 8 line  at 
3.  8 M m for  which  McClatchey  and  Selby  predict  a sea  level  absorption 
coefficient  of  7.3x10  3 km'  ‘ for  a mid-latitude  summer  day  (R.H.  = 80%). 
This  is  shown  as  the  circled  point  on  Fig.  3.  The  dominant  molecular 
absorber  for  this  case  was  found  to  be  HDO. 

Recently  Burch18  has  found  evidence  for  a water  vapor  continuum 
absorption  band  in  the  4 Mm  region.  The  continuum  appears  to  be  due  to 
a pressure  broadening  phenomenon,  with  N2  found  to  be  1 /8  as  efficient 
as  H^O  as  a collision  partner. 

There  is  still  some  uncertainty  in  the  magnitude  of  this  continuum 
absorption  in  that  the  measurements  were  made  at  temperatures  > 338°  K 
and  extrapolated  to  room  temperature.  Nonetheless,  the  water  vapor 
continuum  would  appear  to  be  a prominent  absorption  mechanism  in  the 
4 an  region.  A prediction  of  the  absorption  coefficient  due  to  this  water 
vapor  continuum  is  shown  vs.  R.H.  in  Fig.  3. 

As  can  be  seen  by  reference  to  Fig.  3 molecular  absorption  will 
dominate  aerosol  absorption  for  the  "standard"  continental  aerosol  distribu. 
tion.  On  the  other  hand  much  larger  aerosol  mass  loadings  are  commonly 
observed  in  a variety  of  scenarios,  and  these  in  turn  will  provide  the 


-10- 


om.nant  absorption  loss  mechanism  for  DF  laser  radiation.  This  absorp- 
tion can  be  particularly  severe  under  conditions  corresponding  to  mist, 
heavy  maritime  winds,  etc. 

It  should  be  pointed  out  that  the  main  interest  in  this  phenomenon  is 
not  beam  attenuation,  indeed  scatterinR  from  aerosols  will  attenuate  the 
laser  beam  more  rapidly  than  absorption,  but  rather  the  resulting  effects 
on  beam  quality  due  to  thermal  blooming.  Unlike  molecular  absorption, 
absorption  by  aerosols  results  in  nonuniform  atmospheric  heating  with 
strongly  localized  temperature  gradients.  The  resulting  effects  on  beam 
quality  can  be  significantly  different  from  those  arising  from  uniform  hea‘mg. 
In  the  next  section  a model  is  presented  for  the  prediction  of  the  atmospheric 
perturbations  resulting  from  the  absorption  of  laser  radiation  by  aerosols. 
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III.  MODELING 


The  phenomenon  to  be  modeled  here  is  the  time  and  spatial 
variation  of  the  temperature  and  vapor  fields  about  a water  droplet 
undergoing  laser  irradiation.  This  problem,  in  its  aggregate  form, 
has  been  considered  by  Sutton19  and  Clickler?'°(among  others)  in 
their  analysis  of  laser  propagation  through  fogs  and  clouds,  respectively. 

I he  steady  state  solutions  for  an  individual  droplet,  developed  under  the 
assumption  of  constant  density,  have  been  discussed  by  Williams.  21 
I he  fully  time  dependent  solutions  for  an  individual  droplet  have  not 
appeared  in  the  open  literature. 

I he  phenomenology  of  the  present  problem  is  that  the  aerosol 
droplet  absorbs  energy  and  is  heated.  As  the  droplet  temperature  rises 
above  ambient  the  droplet  conducts  heat  energy  to  the  surrounding  air  and 
undergoes  surface  vaporization.  The  vapor  leaves  the  droplet  with  a 
directed  velocity  and  a p-essure  wave  is  induced  in  the  surrounding  g^s 
in  response  to  the  mass  and  heat  addition.  (If  the  rate  of  energy  absorption 
were  sufficiently  large  the  droplet  would  "explode  ",  creating  a shock  wave; 
the  laser  ^intensity/ particle  size  range  at  winch  this  occurs  has  been  discussed 
by  Sutton  . T he  present  analysis  will  be  limited  to  laser  intensities  and 
droplet  sizes  such  that  the  directed  velocity  of  the  vapor  leaving  the  droplet 
surface  is  much  less  than  the  local  sound  speed),.  An  additional  feature 
which  must  be  included  in  the  analysis  is  the  fact  that  the  droplet  radius 
decreases  with  time  due  to  mass  loss  resulting  from  vaporization. 

‘\ . Gas  Equations 

A number  of  approximations  have  been  made  to  simplify  the  solution 
of  this  problem.  The  most  prominent  of  these  is  that  the  gas  is  taken  to  be 
at  constant  pressure.  As  pointed  out  above,  pressure  gradients  are  induced 
in  the  gas  because  of  the  mass  and  heat  addition  from  the  droplet.  However, 
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if  the  pressure  equilibrates  on  a time  seale  fast  relative  to  tlu.-  times  for 
conduction  and  diffusion  then  the  heat  and  mass  transfer  may  be  assumed 
to  be  occurring  in  a constant  pressure  gas.  The  characteristic  time  for 
" ater  vapor  to  diffuse  a distance  L in  air  is  approximately 

T0  ~ l2/D  (6) 

where  D is  the  diffusion  coefficient.  The  characteristic  time  for  heat 
conduction  in  air  is 

Tc~L2pCp  /k  (7) 

A 

where  p is  the  density,  Cp  the  specific  beat  and  k.  the  thermal  con- 

A A 

ductivity  of  air.  Lastly,  the  characteristic  time  for  pressure  equilibration 
is 

Tp  ~ L/c  (8) 

where  c is  the  sound  speed  in  air.  Evaluating  these  expressions  with 
D = 0.  24  cm  /sec,  PA  = 1.  2 x 10"  * gms/cc,  Cp^  0.  24  cal.  /gm-°K, 

^*2  x 10  cal,  / cm-sec-  K and  c B9  3x10^  cm/sec  results  in 


for  L > i pm.  The  characteristic  droplet  sizes  of  interest  in  this  work  are 
a few  tenths  pm  to  50  Pm,  and  for  the  typical  atmospheric  aerosol  distribu- 
tions of  interest  the  average  spacing  between  droplets,  which  would  be  the 
heating  range  of  interest,  is  approximately  one  mm.  Thus,  it  would  appear 
that  the  approxima  ion  of  a constant  pressure  gas  is  reasonable. 

Another  simplifying  approximation  used  is  that  the  gas  properties 
are  taken  to  be  those  of  clean  air.  Although  this  approximation  would  suggest  that 
the  defining  equations  would  not  be  valid  for  high  water  vapor  concentrations 
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tms  is,  not  the  case,  as  will  be  demonstrated  later  in  this  section.  A’ so, 
the  coefficient  of  thermal  conductivity,  kA  and  the  product  of  the  coefficient 
of  diffusion  and  number  density,  D N,  were  taken  to  be  constant.  These  quan- 
titles  vary  by  < 15%  over  the  temperature  range  of  interest,  0 - 100°  C. 

The  last  simplifications  concern  the  fluid  mechanics.  The  processes 
considered  in  the  analysis  are  thermal  conduction,  diffusion  and  c onvection. 
The  phenomena  of  thermal  diffusion  (Soret  effect)  and  thermo -diffusion 
(Dufour  effect)  have  not  been  included.  Generally,  these  processes  are  in- 
significant in  flow  fields  of  the  type  under  consideration.  Furthermore, 
consideration  must  be  given  to  gravitational  effects.  As  the  gas  surrounding 

the  droplet  is  heated  it  will  rise.  The  characteristic  time  for  the  gas  to 
rise  a distance  L is: 


2L/g 


[^p/p  j ^ 


1/2 


(10) 


*here  g is  the  acceleration  due  to  gravity,  980  cm/sec2,  and  Lp/p  is  the 
fractional  difference  in  density  between  the  heated  gas  and  the  surrounding 
medium.  A strong  lower  bound  on  T ^ corresponds  to  Ap/p  ~ o.  3.  Even  in 
this  limit  it  can  be  show  n that 


td  ~ Tc  <<  Tg 

for  L less  than  on,-  mm.  Of  course  this  effect  might  have  to  be  included  for 
long  time  (t  * 0.01  sec)  irradiation.  However,  in  this  case  it  would  perhaps 

be  more  appropriate  to  treat  the  aerosols  as  an  aggregate  rather  than 
individually. 

Within  the  framework  of  the  above-mentioned  approximations,  the 
general  relationships  describing  the  system  a^e 

a)  Conservation  of  mass 


dp/bt-  -V*  (PV) 
where  V is  the  gas  velocity. 


(ID 
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b)  Euler's  equation 


iV/it  = - V • V V 


c)  Conservation  of  energy 

<J(l/2p  V2  + P Cv  T)/o  t = - V • jp  V (1/2  V2  + CpT)  - kA  V t"|  (13) 

v,here  is  the  specific  heat  at  constant  volume,  and 

d)  The  diffusion  equation 

-*  -»  -* 

0 PV/  dt  V • (D  V Pv  V (14) 

where  Pv  is  the  water  vapor  density  and  D is  the  H^O-air  diffusion  coefficient. 

If  the  droplet  is  taken  to  be  heated  uniformly,  the  system  will  be 
spherically  symmetric  and  it  can  be  shown  that  the  ordered  kinetic  energy 
terms  drop  out  of  the  equations  under  the  condition  of  constant  pressure. 

The  relevant  gas  equations  for  r ■>  R,  where  R is  the  droplet  radius,  are 


k k A T d 

OT/dt  - - V 0T/br  +— — (r2  dT/dr), 

2 or 

r Cp  PM. 

A A 


ax  /6t  = - VoX  /dr  4 k -T-  IP N)  (r2  d X /dr), 
v v 2 o r v 

P r 


dV/dr  - 2 V/r  4 


k kA  A 9 . 2.  . 

1 TT(r  °T/iir)’ 

r cpapma 


where  k is  Boltzmann's  constant,  A is  Avogadro's  number,  P is  the  total 
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pressure  taken  as  one  atm.,  MA  is  the  molecular  weight  of  air,  N is 
the  total  number  density , and  Xv  is  the  mole  fraction  of  water  vapor. 

These  equations  have  been  formulated  in  terms  of  number  density 
rather  than  mass  density,  since  in  a constant  pressure  flow  of  two  Kasos 
with  different  molecular  weights  the  quantity  NT,  rather  than  p T,  remains 
constant.  As  mentioned  earlier,  these  equations  were  derived  under  the 

as  sumption  Xv  1.  In  reality,  the  quantity  Cp  MA  should  be  replaced 
bv  the  expression 

CPA  MA  * V 4 Cpv  Mv  Xv 

ulure  the  subscript  v refers  to  water  vapor,  and  the  diffusion  coefficient 
should  he  replaced  by  that  for  full  binary  diffusion.  However,  over  the 
temperature  and  water  vapor  ranges  of  interest  in  this  problem  these  latter 
two  quantities  are  reasonably  constant  and  well  represented  by  Cp  M and 

D,  respectively.  Thus  Eqs.  (15)  - (17)  are  approximately  valid  for  water 
\ ipoi  concentrations  approaching  unity. 

Droplet  Equations 

In  the  following  analysis  the  temperature  field  across  the  droplet 
is  taken  to  be  constant.  As  will  be  seen  (Sections  IV  and  V)  this  approxi- 
mation is  reasonable  until  the  droplet  temperature  approaches  the  boiling 
point,  if  the  absorption  is  volumetric.  (As  pointed  out  in  Section  II  the 
absorption  will  not  necessarily  be  volumetric.  Specifically,  as  is  evident 
from  Fig.  2,  there  can  be  significant  laser  attenuation  within  the  droplet  for 
large  si,cs.  Ameliorating  this  effect  somewhat  is  the  fact  that  heat  conduction 
within  the  droplet  will  tend  to  smooth  out  temperature  irregularities 
resulting  from  nonuniform  absorption). 

lh'-‘  tja81c  equation  connecting  the  droplet  and  gas  dynamics  is  con- 
servation of  energy  in  the  droplet. 
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(18) 


4/3  TT  R3  I a = 4/3  TT  R3  p Cp  d T /dt 

u D D 

+ 4 TT  R2  kA  (-bT/br  | R) 

2 M 

(-4TTRAH  — — ( N V)  I 

v A ' Mr 


The  LHS  of  Eq.  (18)  is  the  rate  of  absorption  of  energy,  where  I is  the  laser 
intensity,  and  the  RHS  represents  the  energy  deposition  into  droplet  heating, 
conduction  and  vaporization  respectively.  The  quantity  A Hy  is  the  change  in 
enthalpy  required  to  proceed  from  liquid  to  gaseous  state  and  the  subscript  D 
refers  to  droplet  properties.  Note  that  the  temperature  is  continu*  as  across 
the  droplet  surface;  i.  e.  T^  = T (R).  The  terms  involving  the  ordered 
kinetic  energy  of  the  vapor  are  negligible  for  the  present  considerations  and 
have  been  left  out  of  Eq.  (18). 

The  vapor  concentration  at  the  surface  of  the  droplet  has  been 
related  to  the  droplet  temperature  by  the  Clausius -Clapeyron  relationship. 


d In  X 
d~T~ 


P 


AH 


R T 

v 


(19) 


where  R is  the  ideal  gas  constant  per  unit  mass  of  water  vapor.  It  has 

21 

been  demonstrated  that  this  relationship  properly  describes  the  surface 
vapor  pressure  of  a water  droplet  when  the  vapor  velocity  is  less  than  the 
sound  speed  and  the  droplet  radius  is  sufficiently  large  (>  0.05  |jm)  so  that 
surface  tension  effects  may  be  neglected.  In  the  present  analysis  Eq.  (19) 
has  been  used  in  its  integral  form 


Xv  (R) 


(20) 
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where  the  subscript  w refers  to  ambient  conditions.  Eq.  (20)  is  derived  under 
the  assumption  that  A Hv  is  a constant.  Indeed  A Hv  varies  by  less  than  10% 
between  0 - 100CC  and  has  been  taken  to  have  the  constant  value  of  585  cal/gm. 

The  last  equation  connecting  the  droplet  and  gas  behavior  is  the  con- 
servation of  mass  of  the  vapor.  If  a spherical  shell  of  thickness  AR  is 
constructed  about  the  droplet,  then  it  is  required  that  the  rate  of  mass  loss  of 
the  droplet  be  balanced  by  the  rate  of  increase  of  vapor  within  the  shell  plus  the 
rate  of  vapor  leaving  the  exterior  of  the  shell.  With  some  straightforward 
mathematical  manipulation  this  leads  to  the  relationship 

-1 

(NV)  = - (1  - X (R)  ) DN  5X  /Sr  . (21) 

R V V R 

Note  that  Xv  as  given  by  Eq.  (20)  could  exceed  unity  for  sufficiently  high 
temperature.  However,  if  this  occurred,  the  vaporization  rate  as  given  by 
Eq.  (21)  would  oecome  negative.  This  is  a consequence  of  assuming  a constant 
pressure  gas  and  limits  the  droplet  temperature  that  can  be  attained  to  100°C. 

It  was  assumed  in  the  derivation  of  Eq.  (21)  that  the  vaporization  rate  at 
the  surface  is  controlled  by  diffusion.  Actually  vaporization  at  the  surface  is  a 
kinetic  process  obeying  Knudsen's  equation  with  diffusion  taking  control  approxi- 
mately a mean  free  path  (~0.  06  at  1 atm)  away  from  the  droplet.  It  can 
readily  be  show'n^  that  if  the  vaporization  coefficient  for  water,  av,  were  unity 
Eq.  (21)  w'ould  adequately  represent  this  phenomenon  for  particle  radii  > 0.  5|Jm. 

At  the  time  that  this  research  was  in  progress  the  value  of  av  = 1 appeared  to  be 
22-24 

preferred  ; however,  very  recently  a measurement  of  av  = 0.  033  has  been 

2 c 

reported^  in  agreement  with  several  earlier  measurements.  If  this  latter  value 
of  O.  033  is  correct  then  Eq.  (21)  w'ould  provide  an  overestimate  of  the  vaporiza- 
tion rate  for  the  droplet  sizes  considered. 

The  computer  predictions  to  be  presented  in  Section  IV  were  made  using 
Eq.  (21).  The  computer  code  can  be  readily  adjusted  to  include  the  more  general 
vaporization  rate  as  is  discussed  in  Section  V.  In  general,  a smaller  value  of 
av  will  provide  more  heo.t  conduction  for  a given  configuration  at  the  expense  of 
vaporization. 
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C.  Effect  of  Varying  Droplet  Size 


So  far  in  the  analysis  the  effect  of  decreasing  droplet  size  has 
not  been  considered.  This  effect  can  be  uncoupled  because  the  rate  of 
change  of  droplet  size  is  small  compared  to  the  rate  of  change  of  the  gas 
parameters  (because  of  the  large  difference  between  liquid  and  gas 
densities).  Of  course,  it  has  been  implicitly  assumed  that  the  time  varying, 
rather  than  initial,  value  of  R will  be  used  in  the  solution  of  the  equations. 

The  governing  equations  have  been  recast  in  terms  of  the  variable 

Z = R (t)/r  (22) 

In  ter. ns  of  this  variable  the  equations  may  be  written  as 
dT/dt  = R"1  (QZ4  - Z dR/dt)  ST/&Z 


k k A 7 a 7 

+ r — i TR  z » x/a  z2 

CPaPMA 


(23) 


a Xv/dt  = R'1  (QZ4  - Z dR/dt)  a X /dZ 


k DN  _ -2  4 .2  2 

+ — — t r z a x /az 

P v 


(24) 


A A -12  2 

ao/az  = - — A r 1 a2  T/az2 
cp  pma 


(25) 


where 


P M 


d R/d  t 


ApD  k T 


Q (z  l) 


(26) 
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and 


Q = V/Z 


(27) 


The  transformed  boundary  conditions  at  r - R (Z  - 1)  are 


la  = Pn  Cp  bT/bt 
U D 


+ 3 R"2  kA5T/dZ 


Z = 1 


Z r 1 


, M AH  P 

+ 3R  -TTT~~  Q (z  = 1]  ' 


(28) 


q (z  = i)  = r"  1 |1-X„  (Z  = 1)|  axv/az 


[l-Xv(Z  = l)J 


-1 


(29) 

Z^l 


and  the  Clausius -Clapeyron  relationship,  Eq.  (20). 

Note  that  the  quantities  I and  a enter  only  as  a product  in  the 
problem,  Eq.  (28).  As  shown  in  Fig.  2 the  effective  value  of  a as 
determined  by  Mie  theory  can  vary  by  close  to  a factor  of  two  for  particles 
of  radius  between  0.  1 and  100  pim.  Thus  the  absorption  efficiency  of  a 
droplet  can  vary  as  the  particle  vaporizes  (and  thus  decreases  in  size). 
Given  the  uncertainty  in  the  imaginary  part  of  the  index  of  refraction  of 
water  and  the  desire  to  present  generalized  solutions  this  effect  was  not 
included  in  the  computations  of  Section  IV. 
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IV.  MODEL  PREDICTIONS 


No  analytic  solution  has  been  found  for  the  system  of  Eqs.  (20), 

(23)  - (29).  A computer  program  has  been  developed  for  their  numerical 
evaluation  which  employs  a fully  implicit  C rank-Nicolson  finite  difference 
scheme.  A complication  in  using  the  Z co-ordinate  is  that  a constant  grid 
size  in  Z space  results  in  a limited  number  of  grid  points  at  large  r,  where 
the  solution  may  be  most  interesting.  This  complication  has  been  circum- 
vented by  modifying  the  finite  difference  scheme  to  allow  for  a change  in 
grid  spacing  at  small  Z.  In  the  present  mode  of  operation  the  grid  spacing 

for  0.  00  ^ Z ^ 0.  05  is  an  order  of  magnitude  finer  than  that  for  0.  05  £ Z 

^ 1.00. 

A matrix  of  computer  solutions  for  the  Eqs.  (20),  (23)  - (29)  has 
been  completed.  Simultaneously  a steady  state  solution  of  the  equations 
has  been  performed  in  order  to  analytically  predict  the  effects  of  varying 
fundamental  parameters  such  as  particle  radius  and  laser  intensity.  The 
latter  analysis  is  useful  not  only  as  a check  for  the  computer  solutions  but 
also  in  that  it  presents  a broad  overview  of  the  laser-aerosol  interaction. 

A-  Steady  State  Solutions 

In  the  limit  where  the  velocity  term  in  Eqs.  (23)  and  (24)  can  be 
neglected,  the  equations  effectively  uncouple,  and  it  can  be  readily  demon- 
strated that  the  steady  state  solutions  for  A T = T - To,  and  AX  = X - X 
scale  inversely  with  r,  i.  e.  , 00 

A T (r)  = (R/r)  A T (R)  (30 

and 

A Xv(r)  = (R/r)  A Xv  (R)  (31 
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Indeed,  even  the  time  dependent  solutions,  which  may  be  found  in  standard 

texts  {for  example.  Ref.  26)  vary  approximately  as  r ' ^ for  r<  2 (k  t/p  C )^2t 

Thus,  in  the  limit  of  low  velocity  the  gradients  appearing  in  the  equation 

for  droplet  energy  conservation,  Eq.  (28),  may  be  evaluated  by  use  of  Eq.  (30) 

and  (31). 

If  this  is  done,  use  of  Eqs.  (20),  (28)  and  (29)  results  in 


d (ATD)/dt 


R pDCPD 


R2  la 

3 ■ kA  L Td 


AH  M 

D N (1  - Xv  (R)  )_1  (XV(R)  - Xv  ) 


(32) 


where 


X (R)  - X 
v v 


rxp  |_rv(atd  it_) 


Although  this  differential  equation  is  separable,  it  affords  no  simple  analytic 

solution.  One  feature  of  interest  is  the  maximum  value  of  A TD  which  is 

reached  when  the  time  derivative  is  zero.  Equation  (32)  in  this  limit  reduces 

to  a transcendental  equation  for  A T_.  which  ma\  be  evaluated  as  a 

2 ^ * max 

function  of  the  parameter  R la.  The  solution  to  this  equation  is  shown  in 
Fig.  4 over  a range  of  5 orders  of  magnitude  in  the  parameter  R2Ia. 

Shown  for  comparison  are  the  predicted  maximum  droplet  temperatures 
for  the  related  cases  of  "heat  conduction  only"  and  "vaporization  only". 

These  calculations  were  performed  for  initial  conditions  of  T = 293°K 
Xv  = 2.  1 x 10'2. 

V 00 

Note  that  for  the  full  range  of  R2  la  shown  the  AT~  predicted 

D-max 

for  the  case  of  heat  conduction  alone  is  significantly  higher  than  that  pre- 
dicted when  both  heat  conduction  and  vaporization  were  included.  On  the 
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ATq  (max) 


CONDUCTION 


ONLY 


AT  ACROSS 


DROPLET 


fj 


ATD(max) 


VAPORIZATION 


ONLY 


ATq  (max) 


CONDUCTION  AND 


VAPORIZATION 


R Ic^CAL/CM-SEC 


Maximum  Increase  in  Droplet  Temperature  vs.  R2la.  Also 
shown  is  the  fraction  of  absorbed  energy  lost  to  heat  conduction 
and  the  predicted  variation  in  temperature  across  the  droplet. 


other  hand,  the  predictions  for  the  case  of  "vaporization  only"  are  relatively 

close  to  the  full  solution.  The  reason  for  this  is  that  the  major  portion  of 

the  absorbed  energy  goes  into  vaporization  rather  than  heat  conduction.  This 

is  shown  explicitly  in  Fig.  4 where  the  fraction  of  the  absorbed  energy  which 

2 2 

is  lost  to  heat  conduction  is  plotted  vs.  R la.  For  low  R lathis  fraction 

2 

asymptotes  to  ~ 0.  3,  and  falls  off  rapidly  with  increasing  R la  This  rapid 

fall-off  occurs  as  the  droplet  temperature  approaches  the  boiling  point. 

This  is  because  the  vaporization  rate  increases  exponentially  in  this  region, 

whereas  the  rate  of  heat  conduction  approaches  a constant.  The  relative 

value  of  this  fraction  is  determined  by  the  ratio  of  the  last  two  terms  on  the 

RHS  of  Eq.  (32).  For  small  AT^  this  ratio  reduces  to 

D-max 


Power  into  heat  conduction 
Power  into  vaporization 


kARvT-2  A(1  - Xv.» 
*v.  DN(A  Hv)2  Mv 


**  0.  43 


(33) 


in  agreement  with  the  asymptotic  value  shown  in  Fig.  4 (i.  e. , 0.  43/1. 43  «*0.  3). 

Note  that  Fig.  4 is  somewhat  misleading  at  large  values  of  R2  la  in 
that  convective  heat  transfer  will  dominate  conductive  heat  transfer  in  that 
range,  i.  e.  the  water  vapor  is  heated  prior  to  leaving  the  droplet.  The 
fraction  of  the  vaporization  energy  which  goes  to  heating  is 


A T 


h2o 


D 


AHv 


(34) 


which  at  100  C is  Ep  ~ 0.06.  The  fraction  of  absorbed  energy  converted  to 

heat  will  not  fall  below  that  value  given  by  Eq.  (34).  Note  further  that  the 

radial  steady  state  solutions  Eq.  (30)  - (31)  are  not  strictly  valid  in  the 

2 1 

convection  dominated  regime.  Williams  has  discussed  the  steady  state 
solution  for  this  case  under  the  assumption  of  constant  density  and  showed 
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that  the  radial  distribution  is  exponential  with  an  argument  proportional 
to  the  vaporization  rate  over  radius.  Nonetheless  for  the  cases  corresponding 
to  high  vaporization  considered  in  the  present  work  the  computed  racial 
derivatives  of  the  temperature  and  water  vapor  distributions  at  r R were 
typically  within  a factor  of  two  of  those  predicted  from  Eqs.  (30)  (31). 

Thus  it  is  concluded  that  the  present  steady  state  analysis  should  adequately 
describe  the  gross  phenomenology  of  the  convection  dominated  regime. 

One  last  point  with  reference  to  Fig.  4 is  that  although  the  solution 
for  the  case  of  "heat  conduction  only"  is  linearly  proportional  to  laser 
intensity  this  is  not  the  case  for  the  more  general  solution.  As  vaporization 
begins  to  strongly  dominate,  the  temperature  increase  scales  more  weakly 
with  laser  intensity.  This  occurs  approximately  at  an  R2  la  of  1.5  x 10~2 
cal/cm- sec  which  for  a laser  intensity  of  1 04  watts/cm2  and  an  a of  112  cm-l 
corresponds  to  a droplet  size  of  ~2pm. 


The  fact  that  vaporization  dominates  heat  conduction  as  an  energy 
loss  mechanism  allows  an  approximate  analytic  formulation  for  several  of 
the  variables  of  interest  in  the  problem.  In  particular,  in  steady  state 
it  can  be  readily  shown  from  Eq.  (32)  that  the  rate  of  mass  ioss  from  the 
droplet  can  be  expressed  as 


dm/dt  - - 4/3  tt  R3  la /A  H 

v 

or  in  terms  of  the  droplet  radius 


(35) 


dR/dt  = - RI  a/(3  A pD)  (36) 

From  Eq.  (36)  it  can  be  seen  that  the  characteristic  e- folding  time  for 
droplet  radius  is 


TR  = 3 AHvVIa 


(37) 
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which  is  inversely  proportional  to  la  and  independent  of  initial  particle  si/e. 
Lastly,  the  gas  velocity  leaving  the  droplet  surface  is  given  by 


V (R)  60  MA  R la/  (3MvPi  Hy)  (38) 

2 

Note  this  latter  quantity  scales  as  R la  rather  than  R la,  and  thus  is  less 
sensitive  to  changes  in  radius  than  the  droplet  temperature. 

The  relationships  (33)  - (38)  are  all  valid  after  the  droplet  nas  reached 

its  steady  state  temperature.  It  can  be  shown  that  prior  to  that  time  the 

major  portion  of  the  absorbed  energy  is  consumed  in  heating  the  droplet.  It 

was  found  through  computer  solutions  of  Eq.  (32)  that  the  o- folding  time  for 

the  auantity  A - A was  given  bv 

D-max  D 


0.  65  AT 


td 


D-max  L 


la 


(39) 


to  within  20%  for  R^  la  between  10 

not  independent  of  particle  size  since 
2 -2 

for  R la  < 1.5x10  cal/cm-sec 


- 10  cal. /cm- sec.  Note  that 

2 

A T_  is  a function  of  R lo.. 

D-rnax 


is 

Indeed 


4Tn  ~ 1.  6 x 103  RZ  la  (40) 

D-max  ' ' 

whence 

TD~  103  PL  CpL  R ; R Id  < 1.5x10  cal/ cm- 3ec  (41) 

independent  of  la.  As  a numerical  example,  for  a droplet  of  radius  1pm  the 

_ 5 

characteristic  time  to  steady  state  as  determined  from  Eq.  (41)  is  10  seconds. 

For  values  of  R la  larger  than  1.  5 x 10  cal/cm-sec  the  characteristic 

heating  time  must  be  determined  from  Eq.  (39)  in  conjunction  with  the  values 

of  A T_.  shown  in  Fig.  4. 

D - max 
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A last  point  before  leaving  the  discussion  on  steady  state  solutions 
concerns  the  importance  of  temperature  gradients  within  the  droplet.  If 
one  considers  volume  absorption  followed  by  heat  conduction  as  the  only 

phenomenon  occurring  within  the  droplet,  the  defining  steady  state  energy 
equation  becomes 

4/3  TT  r^  la  = - 4 TT  r b T/br  (r  < R)  (42) 

where  kD  is  the  coefficient  of  thermal  conductivity  of  liquid  water,  taken  as 
1.4  x 10  cal /cm- sec -°K.  Eq.  (42)  requires  that 


^ ^ R la 

C - TS  = 6 k (43) 

where  the  subscripts  C and  S refer  to  conditions  at  the  droplet  center  and 
surface,  respectively.  This  temperature  variation  is  also  shown  in  Fig.  4 
and  can  be  seen  to  be  significantly  lower  than  &TD_max  until  the  droplet 
approach-  the  boiling  point.  Of  course,  Eq.  (43)  is  not  necessarily  realistic 
in  that  region  since  it  requires  droplet  superheating.  In  reality,  vaporization 
could  occur  within  the  droplet  and  original  nuclei  effects,  etc.  could  be  quite 
important.  These  points  are  discussed  in  Section  V.  In  any  event,  it  is 
clear  that  the  approximation  of  a constant  temperature  field  across  the 
droplet  is  not  appropriate  at  large  values  of  R2  la. 

Time  Dependent  Solutions 

Computer  calculations  have  been  performed  for  droplet  radii  between 
0.  5 - 50  pm  and  for  values  of  Ia  = 27,  81,  270,  810.  2700  kcal/cc-sec. 

(For  a nominal  value  of  a = 113  cm"1  these  correspond  to  I = 1,  3,  10,  30, 
100  kwatts/cm  ).  Only  a few  predictions  have  been  made  for  the  largest  and 
smallest  droplet  sizes  inasmuch  as  these  require  the  largest  amount  of 
computer  time. 
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Shown  in  Figs.  5-9  are  the  rise  in  droplet  temperature  vs. 

time  for  an  I a of  270  kcal/cc-sec  (I  = 1 04  watts/cm2  for  an  a of 
- 1 

113  cm  ) and  for  particle  radii  of  0.  5,  1.,  3.,  5.  and  10.  Mm.  The 
slope  discontinuities  appearing  on  the  curves  are  not  real  but  are  a 
result  of  the  plotting  technique.  As  can  be  seen  the  time  to  reach  steady 
state  increases  with  increasing  particle  size  as  expected  from  Eq.  (39)  and 
(41).  The  decrease  in  droplet  temperature  with  increasing  time  results 
from  the  decrease  in  radius  due  to  vaporization.  The  predicted  droplet 
temperatures  for  times  corresponding  to  the  point  of  maximum  tempera- 

and  beyond  are  in  good  agreement  with  the  steady  state  predictions  shown 
in  Fig.  4. 

The  radial  profiles  of  the  temperature  increase  in  the  gas  for 

these  cases  are  shown  in  Figs.  10  - I 4 at  times  of  1 , 3 and  5 msec. 

These  profiles  scale  inversely  with  radius  out  to  distances  roughly  given 

1/2 

by  r ~ 2 (k  At/p  Cp  ) in  agreement  with  the  analytic  solution  for  the 

A 

case  of  heat  conduction  with  no  vaporization.  (The  magnitude  of  the  increase 
is  of  course  lower  than  the  heat  conduction  case  as  is  evident  from  Fig.  4). 
The  decrease  in  temperature  at  small  radius  with  increasing  time  is 
caused  by  the  decreasing  droplet  temperature.  At  the  larger  distances 
the  temperature  is  highest  at  the  longest  time  because  of  the  finite  time 
required  for  heat  conduction. 

The  radial  profiles  of  the  increase  in  water  vapor  mole  fraction 
in  the  gas  are  quite  similar  in  form  to  those  for  the  temperature  increase. 
Examples  of  this  are  shown  in  Figs.  15  and  16  for  particle  radii  of  1 and  10 
^ m and  an  I Ct  of  270  kcal/ cc - s ec.  The  profiles  scale  inversely  with 
radius  out  to  approximately  the  same  distance  as  the  respective  profiles 
for  the  temperature  increase  in  the  gas  as  shown  in  Figs.  11  and  14.  This 
is  because  the  characteristic  time  for  diffusion  is  approximately  the  same 
as  that  for  heat  conduction  (See  Eqs.  (6)  and  (7)  ). 
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Similar  predictions  for  an  i « of  2700  kcal/cc-sec  (I  - 10^  watts/cm 
for  an  a of  i 1 3 cm  ) and  for  particle  radii  of  0.  5,  l.,3.  and  5.  Mm 
are  shown  in  Figs.  17-24.  The  rise  in  droplet  temperature  vs.  time  for 
various  droplet  sizes  is  shown  for  this  I a in  Figs.  17-20.  Note  the 
time  scale  to  reach  steady  state  again  increases  with  radius.  At  the 
larger  droplet  sizes  these  times  are  considerably  faster  than  those 
for  the  same  droplet  size  at  an  1 a of  270  kcal/cc-sec.  This  is  as 

expected  from  Eq.  (39).  The  fall-off  in  droplet  temperature  with 

increasing  time  is  more  rapid  in  this  case  than  at  the  lower  la  since 

the  characteristic  time  for  droplet  shrinkage  scales  inversely  with  1 a 
(Eq.  (38)  ). 

The  increase  in  temperature  of  the  gas  vs.  radius  is  shown  in 
Fig.  21-24  at  times  of  0.  2,  0.  6 and  1 . 0 msec.  Departures  from  an 
inverse  radial  dependence  are  obvious  in  these  profiles  particularly  for 
the  larger  droplet  radii.  The  corresponding  predictions  for  the  increase 
in  water  vapor  mole  fraction  exhibit  a similar  behavior. 

Predictions  of  the  cumulative  energy  distribution  vs.  time  are 
shown  in  Figs.  25-29  for  a droplet  radius  of  5 Hm  and  for  five  values 
of  I a lying  between  27  to  2700  kcal/cc-sec.  Shown  are  the  total  absorbed 
energy,  the  amount  of  energy  converted  to  heat  of  the  surrounding  gas  by 
conduction  and  by  convection  and  that  lost  to  uaporization.  Note  that  the 
five  figures  are  not  all  plotted  on  the  same  time  scale. 

The  first  thing  that  can  be  seen  is  that  the  total  energy  absorbed 
does  not  scale  linearly  with  1 a over  the  time  ranges  shown.  The  reason 
for  this  is  that  the  rate  of  droplet  shrinkage  increases  with  increasing 
I a and  thus  the  effective  absorption  coefficient  decreases.  Also  the 
percentage  of  the  total  absorbed  energy  which  is  conducted  away  as  heat 
energy  decreases  with  increasing  la  as  expected  from  the  steady  state 
analysis  discussed  in  Section  A.  Indeed  the  heat  flow  resulting  from  con- 
vection exceeds  that  from  conduction  at  an  I a of  2700  kcal/cc-sec. 
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It  should  be  emphasized  that  the  prediction  that  vaporization, 
rather  than  heat  conduction,  is  the  dominant  energy  loss  mechanism  is 
strongly  dependent  upon  the  assumption  that  the  evaporation  coefficient 
for  H^O  is  unity.  If  it  were  0.  033,  as  mentioned  in  Section  II  conduction 
could  dominate  vaporization.  This  point  is  discussed  in  some  detail  in 
the  next  section. 
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FIGURE  CAPTIONS 


Figs.  5-9  The  Rise  in  Droplet  Temperature  Above 

Ambient  vs.  Time.  lit  2.  7 x 10^  cal/cc-sec, 

R = 0.  5,  1.  , 3.  , 5.  , 10.  pni,  respectively. 

Ambient  conditions  arc  T = 293°K,  Xv  2.  1 x 10-2. 

Figs.  10  - 14  The  Rise  in  Gas  Temperature  Above  Ambient  vs. 

Radius  at  Times  of  1,  3,  5 msec.  Initial  conditions 
the  same  as  in  Figs.  5-9,  respectively. 

Figs.  15  & 16  The  Rise  in  Water  Vapor  Mole  Fraction  Above 
Ambient  vs.  Radius  at  Times  of  1,  3,  5 msec. 

Ia  - 2. 7 x 10^  cal/cc-sec,  R = 1.,  10.  p m 
respectively.  Ambient  conditions  are  T - 293°K, 
Xv  = 2.  1 x 10"2. 


Figs.  17-20  The  Rise  in  Droplet  Temperature  Above  Ambient 
vs.  Time.  Ia  = 2.  7 x 10^  cal/cc-sec,  R = 0.  5, 

1.  , 3.  , 5.  MTn  respectively.  Ambient  conditions  are 
T 293°K,  Xv  r.  2.  1 x 10'2. 

Figs.  21-24  The  Rise  in  Gas  Temperature  Above  Ambient 

vs.  Radius.  Initial  conditions  are  the  same  as  in 
Figs.  17  - 20,  respectively. 

Figs.  25  - 29  Cumulative  Energy  Deposition  vs.  Time. 

R - 5 pm,  Ia  = 2.  7 x 104,  8.  1 x 104, 

2.  7 x 105,  8.  1 x 105,  2.  7 x 10<>  cal/cc-sec, 
respectively.  Shown  are  the  total  absorbed 
energy  and  the  amounts  lost  to  heat  conduction, 
heat  convection  and  vaporization. 


-32- 


Droplet  Temperature  v Time 
IflLPHA  - 270  kcal/cc-sec  , R - 0.5  micron 
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Droplet  Temperature  v Time 
I ALPHA  - 270  kcal/cc-sec  , R - 3 microns 


Droplet  Temperature  v 
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Fig.  13 
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Vopor  Concentration  Increase  v Radius 
I ALPHA  - 270  kcal/cc-sec  , R - 1 micron 


.j 

Droplet  Temperature  v Time 
I ALPHA  “ 2700  kcal/cc_sec  , R - 1 micron 


Fig.  13 
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Temperature  Rise  v Radius 


I ALPHA  - 2700  kcat/cc~sec  , R - 1 micron 
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V.  "NON- IDEAL"  DROPLET  EFFECTS 


The  modeling  discussed  in  the  previous  sections  has  been  limited 
to  the  ideal  case  of  a pure  water  droplet  undergoing  volume  absorption 
As  alluded  to  in  Section  II  there  are  a number  of  effects  pertaining  to 
real  aerosol  properties  which  either  were  not  included  in  the  model  or 
were  tncluded  in  an  approximate  manner.  Some  of  these  effects  are 
considered  in  detail  in  this  section.  These  include  consideration  of 

non-untform  absorption  and  heating,  vaporisation,  solution  effects  and 
droplet  ebullition. 


A.  Nonuniform  Heating 
1 . Ray  Optics 

The  calculations  discussed  in  Sections  III  and  IV  have  all  been 
based  on  the  concept  of  'Volumetric  absorption",  which  is  something 
of  a mathematical  fiction.  From  a ray  optics  viewpoint,  it  assumes 
t all  rays  pass  through  the  irradiated  volume  without  deflection 
and  with  negligible  attenuation;  the  energy  deposition  along  a ray 

bundle  of  area  dA  is  lad  A dr,  or  la  watts  per  unit  volume  at 
all  pomts  within  the  irradiated  volume,  regardless  of  its  shape. 

Ray  optics  can  be  used  to  examine  the  water  droplet  more  realistically 
choosing  a droplet  sire  for  which  the  theory  is  applicable,  without  at 
this  point  specifying  how  large  that  size  must  be.  From  the  sketch 
below  it  is  clear  that  most  of  the  rays  have  longer  path  lengths  within  the 
droplet  than  they  would  have  if  they  were  not  deflected  by  refraction 
a.  the  droplet  surface.  For  refractive  index  n,  Snell's  law  yields 
s.n  « = „ sin  <S,J and  if  an  absorption  coefficient  a cm’1  is  assumed 

then  I(s)=Ioe  , where  s is  measured  along  the  ray  from  its 

intersection  with  the  droplet  surface.  The  resultant  energy  deposition 
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rate  within  the  droplet  relative  to  the  volumetric  absorption  rate  can  be 
expressed  by  the  ratio 

>.  R y»2  R cos  i 

J 2 TTy  dy  Jn 


rQ  2 TTy  dy  JQ  Id  S 

R /*R  2 TTy  dy  (I  2R  cos  0) 
JQ  u 


(44) 


Making  the  substitutions  z = y2  / R2  = sin2  q,  and  cos  4>  =>/l  - z/n2 
one  obtains 

A = 0 a L Jo 


fld z 2Ry/l  - z 

Jo 

f0  - exp^-2  aR^l  - z/n2)J 


dz 


2aR y \/l  - z dz 


4aR  1 " 2 ( 2aR  ) K1  +Z0Ryfl  - l/n2)exp(-2aR>/l  - 1/n2) 

- (1  + 2a R)  exp(-2aR)| 

This  expression  has  several  interesting  limiting  values.  For  instance, 
when  n = 1, 


(45) 
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3 

4aR 


(1+20R)  exp(-20R) 
' 2a2  R2 


] 


1 - 3aR/4  as  0R->  0 . 


(46) 


Thus  ray  optics  recovers  the  volumetric  rate  (A  = 1)  when  a R -»  0, 

R 

and  the  correction  for  finite  absorption  is  obtained  in  the  case 
where  rays  are  not  deflected. 

For  finite  n,  and  oR-^0,  Eq.  45  takes  the  form 

AR  = "2  [*  - (l  - l/n2)3/2]  , (47, 

yielding  AR  = 1.  275  when  n = 1.  364  . Finally,  for  large  values  of 
OR,  AR-»3/4dR,  which  is  independent  of  n.  This  implies  a 
correction  to  the  volumetric  rate;  i.  e. 

4 3 o 

total  energy  deposited  = — TTR  I ctA  = TTR  i . 140, 

->  o R o ' ' 

Thus  for  large  a the  incident  radiation  is  entirely  absorbed  within  the 
droplet  surface,  and  refraction  is  irrelevant. 

For  n=  1.  364  and  a = 1 12.  4 cm"1,  Eq.  45  has  been  used  to 
evaluate  AR  as  a function  of  R,  and  values  of  A a are  shown 
in  Fig.  2.  ARa  is  thus  an  effective  absorption  coefficient  and  it  may  be 
compared  with  the  absorption  coefficient  obtained  from  Mie  theory, 
as  discussed  in  Section  II,  and  shown  in  Fig.  2.  It  is  clear  that  the  ray 
optics  value  of  effective  absorption  coefficient  is  excellent  for  R>50  Mm, 
but  rather  seriously  in  error  (as  would  be  expected)  for  R < 20  Mm. 


Knowing  the  focusing  effect  of  a spherical  refracting  surface,  one 
can  be  sure  that  the  heating  within  a droplet  of  R > 50  is  highly 
nonuniform.  For  smaller  droplets,  Mie  theory  establishes  that  the 
total  absorption  becomes  much  higher  than  the  volumetric  rate,  but  no 
statement  can  be  made  about  the  spatial  nonuniformity  of  the  energy 
deposition  at  this  time.  This  can  only  be  established  by  using  Mie 
theory  to  evaluate  the  fields  and  their  derivatives  within  the  droplet, 
and  this  effort  does  not  appear  justified  at  present.  However,  it  is 
worthwhile  to  revert  to  the  ray  optics  and  ascertain  the  local  energy 
deposition  rate  as  a function  of  position  within  the  droplet.  Establishing 
the  location  of  contours  of  constant  absorption  rate  is  a straight-forward 
bookkeeping  problem  on  a computer,  and  two  cases  have  been  evaluated. 
Figure  30  shows  the  location  of  contours  for  M = 3,  4,  5 and  10,  where 
M is  the  ratio  of  local  absorption  rate  to  the  volumetric  absorption  rate. 

Along  the  droplet  centerline  M increases  monotonically  from  1.  0 

near  the  front  face  to  slightly  over  4 at  the  rear  face.  Higher 

absorption  rates  are  found  in  a cusp-shaned  region  where  ray-crossing  occurs 

Beyond  this  region  lies  almost  20%  of  the  sphere  volume  which  receives 

no  energy  input  at  all.  The  results  are  displayed  more  quantitatively 

in  Fig.  31.  The  curve  of  M versus  p shows  what  percentage  (P) 

of  the  sphere  volume  is  contained  within  any  given  M-value  contour. 

Corresponding  to  any  P is  a value  N which  represents  the  mean  energy 

deposition  rate  (relative  to  the  volumetric  rate)  for  the  entire  volume 

within  the  M(P)  contour.  Thus  the  curves  show  that  2%  of  the  volume 

lies  within  the  M = 5 contour,  while  the  mean  rate  within  this  2% 

volume  is  N = 8.  8 times  volumetric. 

As  long  as  clR  , 01  there  is  negligible  attenuation  along  a ray, 
and  the  results  labeled  "no  absorption"  are  appH cable  for  all  values  of 
R for  which  ray  optics  is  valid.  Unfortunately,  for  R>50  [Jm  the 
attenuation  is  appreciable,  so  a separate  calculation  has  been  carried 
out  for  the  R = 50  H m case,  and  plotted  in  Fig.  31.  although 
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the  focusing  effects  are  less  dramatic,  one  still  finds  1%  of  the 
sphere  volume  within  the  M = 3.  5 contour,  and  an  average  rate  of  5 
times  volumetric  within  that  contour. 
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2*  Temperature  Gradients  Within  a DrnpW 

Although  at  this  point  no  sin  pie  way  of  determining  the  magnitude 
of  the  nonuniformities  in  heating  rates  within  the  smaller  droplets 
(R<20  urn)  has  been  established,  it  would  appear  that  nonuniformities 
Wtll  exist.  Thus  it  is  important  to  assess  the  effects  of  internal 
temperature  gradients  caused  by  these  nonuniformities. 

It  has  already  been  pointed  out  that  even  the  uniform  absorption  requires 
a nonuniform  temperature  distribution  within  the  droplet,  since  heat 
must  flow  toward  the  droplet  surface  for  removal  by  conduction  and 
vaporization.  For  a heating  rate  of  1q  ^ the  steady-state  difference 

in  temperature  between  droplet  center  and  surface  was  shown  to  be 
I a k 


ATo  = ' »hl!re  kD  "1-4  x 10'3  cal/cm-sec°  K is  the 

thermal  conductivity  of  water.  The  value  of  AT  is  plotted 

in  Fig.  4 as  a function  of  1q  OR2.  In  the  ^ examples  which  foll<w 
the  effect  of  nonuniform  heating  on  AT  = T _ T 

examined,  and  the  results  expressed  at  A t/At"  For^mplicity, 

both  cases  considered  preserve  the  spherical  symmetry,  and  steady- 
state  conditions  are  assumed. 

For  the  first  example  a uniform  heating  rate  la  cal/cm3  - sec 
Within  a radius  r < R is  assumed,  while  the  remainder  of  the  droplet 

is  taken  to  be  heated  at  the  primary  rate  1^.  Then  it  can  be  shown 


(49) 


-f»'1  - 


and  this  function  is  plotted  in  Fig.  32.  The  curves  indicate  that  if 

~ 1%  of  the  central  volume  of  the  droplet  had  a heating  rate  of  10  times 

the  volumetric  rate,  then  the  AT  value  would  be  doubled.  This  localized 

heating  need  not  be  restricted  to  the  focusing,  of  course;  this 

analysis  is  equally  valid  for  the  case  of  an  insoluble  nucleus  having  a 

large  absorption  coefficient  within  the  droplet  for  which  a » a . 

Fig.  30  showed  that  the  focusing  effects  generated  high  heating  rates 

much  closer  to  the  surface  than  to  the  center.  Search  for  a solution  of 

an  asymmetric  heating  problem  of  this  type  is  not  justified  at  this 

time,  but  a spherically  symmetric  case  of  some  relevance  can  be  readily 

set  up.  This  is  illustrated  in  Fig.  33,  where  now  the  higher  heating 

rate  (la)  is  confined  to  a spherical  shell  r^  <r  < r^,  while  the 

basic  rate  I a applies  to  the  regions  0<  r<  r and  r < r<  R. 
o o 12 

In  this  case  the  steady  state  solution  is 


at 

AT 


la 

I a 
o o 


(50) 


The  effect  of  varying  the  position  of  the  r^r^  shell  is  best  displayed 

by  keeping  the  shell  volume  constant.  Thus  Fig.  33  shows  values  of 

AT/AT  plotted  versus  r,  /R.,  for  shell  volumes  of  5,  10  and  15% 
o 1 

of  the  droplet  volume.  As  r,  ^0  the  values  of  AT/AT  are  identical 

1 3 

with  those  obtained  from  Fig.  32  with  (r/R)  = .05,  . 10  and  . 15 

respectively.  Thus  the  curves  of  Fig.  33  show  the  decrease  obtained 

in  AT /AT  as  the  region  with  higher  heating  rate  is  displaced  from 
o 

the  center.  Nevertheless  it  is  clear  that  the  possible  increases  in 

T - T are  still  appreciable,  and  this  can  be  an 

center  surface 

important  factor  once  T J approaches  the  boiling  point,  as  will 
r center 

be  discussed  later. 
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Differential  Heating  Within  Spherical  Shell.  The  volumetric 
absorption  rate  la  applies  to  the  region  rj<r  r?.  Values  of 
LY  Tr  - Tc  are  shown  as  functions  of  rj  for  shell  volumes  of 
5,  10  anti  15%  of  the  total  sphere  volume. 


B.  Vaporization  Coefficient 

In  all  the  cases  shown  in  this  report  the  mass  loss  calculations 
are  based  on  diffusion  with  a driving  vapor  pressure  concentration  at 
the  droplet  surface  assumed  to  be  the  saturated  vapor  pressure  at  the 
surface  temperature.  This  assumption  has.  limitations  which  have  been 
known  for  many  years.  To  assess  the  effect  of  these  limitations,  one 
may  examine  the  steady  state  diffusive  evaporation  of  a droplet, 
following  the  analysis  given  by  Fuchs.  The  straightforward  Maxwell 
solution  for  vaporization  in  the  absence  of  temperature  gradients  gives 
the  diffusive  flux  as 


J - 4 tt  r D | c (r)  - c^  | g/sec  (51) 

2 

where  D is  the  diffusion  coefficient  (cm  /sec),  and  water  vapor  concen- 

3 

trations  c (r)  and  c^  are  in  g/cm  . Near  the  droplet  the  flux  may  also 
be  written  as 


J 


4 rr  R 


(co  ‘ 


Cj)  va 


(52) 


where  c^  is  the  concentration  at  the  surface  R,  c^  the  concentration  at  a 

distance  A above  the  surface,  V = VkT/2H  m cm/sec,  and  a is  the 

v 

vaporization  coefficient  for  the  liquid,  usually  set  equal  to  the  condensation 
coefficient  for  the  vapor.  In  the  region  R < r < R + A,  where  A is  roughly 
two  or  three  times  the  gaseous  mean  free  path,  the  concentration  gradient 
is  much  steeper  than  in  the  region  where  steady  state  diffusion  is  established. 
Since  Eq.  (51)  gives  J = 4TT  (R  + A ) D (c^  - c^  ),  we  may  use  Eq.  (52)  to 
eliminate  c^  and  obtain 


J 


4 TT  R D (c  - c ) 
o 00 


(53) 
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Thus,  the  Maxwell  flux  obtained  by  Betting  r = R in  Eq.  (51)  is  reduced  by 
the  factor  in  the  denominator  of  Eq.  (53).  Now  at  sea  level  the  mean  free 
path  is  - 0.06  pm.  so  that  the  term  R/(R  + A ) does  not  differ  appreciably 
from  unity  unless  R < 0.5(1  m,  which  is  approximately  the  lower  limit  of 
interest  for  the  present  study.  The  other  term  may  be  approximated  as 


D 


Rva 


6 a 


(104/R) 


(54) 


and  at  the  start  of  the  present  study  it  was  felt  that  this  term  could  be 


eglected  for  R > 0.  j (-i  m,  since  the  most  recent  experimental  evidence22"24 


indicated  that  a ^ 1.  Under  these  circumstances  the  two  terms  in  the 


denominator  of  Eq.  (53)  tend  to  be  s elf- compensating,  and  their  sum  remains 
close  to  unity. 


However,  recent  work25  has  tended  to  revive  the  justification  for 


using  av  ~ . 033  for  water  condensation,  so  that  the  term 


D 


R va 


5 (10"4/R) 


(55) 


which  does  not  become  negligible  until  R > > 5 pm.  Even  for  much  larger 

drops,  when  the  reduction  factor  approaches  unity,  one  still  finds  a signi- 
ficant effect  on  the  vaporization  process  if  a < < 1. 


As  pointed  out  in  Section  IV  the  diffusion  solution  is  only  valid  for 
values  of  R la  where  vaporization  is  not  strongly  dominant  and  convection 


terms  may  be  neglected.  For  small  droplets  which  do  not  reach  the  boiling 


point  (as  an  example  from  Fig.  4 for  droplets  with  R < 3 Mm  for  a nominal 
value  of  la  = 2.  7 x 105  cal/cm3-sec)  the  diffusion  solution,  Eq.  (53),  is 


appropriate.  From  Eq.  (55)  it  can  be  seen  that  the  vaporization 
rate  will  be  considerably  smaller  than  that  predicted  in  Section  IV  if 
a = 0.  033  rather  than  unity.  Indeed  the  quantity  1 + D/(R  v a ) may  be 


envisioned  as  a correction  to  the  diffusion  coefficient  at  the  droplet  surface,  i.  e. 


-b'l- 


■ ' 


• 


R v a 

V 

In  Section  IV  it  was  shown  that  the  ratio  of  the  absorbed  energy  going 
into  conduction  to  that  going  into  vaporization  was  inversely  proportional 
to  the  diffusion  coefficient  at  the  surface  for  droplet  sizes  with  steady 
state  temperatures  < 40°  C (see  Eq.  (33)  ).  If  Eq.  (56)  is  used  for  this 

diffusion  coefficient,  witho^  = 0.033,  Eq.  (33)  would  predict  that  con- 
duction strongly  dominates  vaporization  for  these  particle  sizes.  As 

discussed  in  Section  IV  the  opposite  would  be  true  if  a were  unity.  A 

corollary  to  this  increased  conduction  would,  of  corn-  e,  be  that  the 
predicted  steady  state  temperatures  of  these  droplets  would  approach 
those  predicted  for  the  case  of  heat  conduction  alone  as  shown  in  Fig.  4. 

An  additional  consequence  of  Eq.  (52)  occurs  at  large  droplet  sizes 
or  laser  intensities  where  vaporization  dominates.  It  was  pointed  out  in 
Section  IV  that  in  this  limit 


J = - dm/dt  = 4/3  tt  r3  ia  /ah  , 

v ' 

Since  it  is  known  from  Eq.  (52)  that  this  mass  efflux  is  kinetically  limited, 
Eq.  (34)  must  put  a constraint  on  the  surface  vapor  concentration.  The 
modeling  presented  so  far  was  developed  under  the  constraint  of  a constant 
pressure  gas.  It  is  a necessary  but  obviously  not  sufficient  criterion  that 

<Co  " T < *0  dynes/cm2  (i 

for  the  gas  to  be  at  constant  pressure.  (Actually  this  quantity  should  be 

considerably  less  than  106  dynes/cm2  since  = R + A ~ R).  If  Eq.  (52) 

and  Eq.  (34)  are  equated  it  can  be  readily  seen  that  the  constraint  (57)  will 
only  be  satisfied  if 
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R la 


3 x 108  a AH 

v 

1 / 7 

(2  TT  R T) 
v 


(57) 


or  for  = 0.033,  and  T = 373°  K, 


R la  < 600  cal/cm2-sec  ^ 

For  a nominal  value  of  la  = 2.  7 x 105  cal/cc-sec  this  corresponds  to 
R < 22  Mm.  Thus,  it  would  appear  that  ifav  = 0.  033,  significant  over- 
pressure (probably  leading  to  shock  formation)  could  occur  at  relatively 
moderate  particle  sizes  and  laser  intensities.  The  isobaric  vaporization 
modeling  developed  in  this  report  would  not  be  valid  in  such  situations. 

C.  Solution  Droplets 

All  of  the  work  presented  above  has  been  explicitly  concerned  with 
pure  H20  droplets,  although  it  is  well  recognized  that  all  naturally  occurring 
aerosol  droplets  contain  original  nucleus  material,  which  may  be  entirely 
in  solution  or  remain  partially  insoluble.  Junge  and  McLaren28  claim  that 
very  few  aerosol  nuclei  contain  less  than  25%  soluble  material,  and  cite 
measurements  for  which  70%  of  the  cases  contained  more  than  50%  soluble 
material. 

1.  Droplet  Growth 

Much  work  on  growth  of  solution  droplets  has  been  carried  out  by 

cloud  physicists,  and  an  extensive  literature  is  largely  applicable  to  the 

present  situation.  The  equilibrium  vapor  pressure  p'  at  the  surface  of  a 

r 

droplet  of  radius  R is  expressed  relative  to  the  vapor  pressure  p over  a 

*00 

plane  H^O  (pure)  surface  at  the  same  temperature.  It  will  be  shown  that 
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(a„.)  exp 
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R P'  R T 
D v 


where  &w  is  the  activity  of  the  solvent  and  o',  p'D  are,  respectively, 
the  surface  tension  and  density  of  the  solution. 

Variants  of  this  equation  exist  in  the  literature:  the  present 
discussion  combines  the  approaches  of  Mason10  and  Low  29.  The 
change  in  free  energy  when  an  elemental  mass  dm  of  H^O  i.  transferred 
from  a solution  droplet  to  the  plane  water  surface  is 


- 0 1 d A - n dv  = dm/p1^  (2a  1 /R  - fl  ) 


where  17 

of  the  activity  of  the  so 

■’'•***  *o  o lUIlLtlUIl  ui 

concentration.  The  vapor  pressure  is  obtained  by  noting  that  the 
relationship 


'Rv  TPw  ln  aw  is  the  osmotic  pressure  expressed  in  terms 
lvent,  which  itself  in  turn  is  a function  of  solute 


£ G = dm  RvT  ln  (pWp J 

expresses  the  change  in  free  energy  to  vaporize  water  at  p'^,  compress 

to  Pco  , and  condense  back  to  liquid.  Equating  the  two  expressions  for  A G 
results  in 


(60) 


(61) 


ln  P'/Pco  = ~T  In 
D 


2 a ' 

a + 

w R P'  R T 


D v 


(62) 


which  is  equivalent  to  Eq.  (59). 

The  solvent  activity  a is  generally  tabulated  as  a function  of 

molality  of  the  solution.  Low  9 provides  tables  for  eight  common  electro- 

lytos  which  are  found  in  aerosols.  To  obtain  the  molality  for  a solution 

containing  m^gni^f  solute  of  molecular  weight  M'  one  first  obtains  the 

H2°  mass’  T R PD  " m'*  The  molality,  ms  (moles  of  solute  per  kg 
°1  > Is  then 
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m 


1000  m1  /M1 


(63) 


s 4 3 

3 n R P'D-  m* 

An  important  aspect  of  the  application  of  Eq.  (59)  is  the  fact 

that  for  a given  mass  of  solute  one  can  in  general  find  a droplet  radius 

consistent  with  a value  of  p'^  which  yields  pWp^  equal  to  the  prevailing 

relative  humidity.  Note  that  this  can  never  happen  for  pure  HO  droplets 

for  which  a =1,  since  an  equilibrium  droplet  then  always  requires 
w 

supersaturated  surroundings.  This  was  to  some  extent  a drawback  in  the 

pure  H^O  calculation  of  drcplet  evaporation;  if  the  surface  tension  term 

is  neglected  and  the  R.H.  maintained  at  100%  a pure  H^O  droplet  remains 

in  equilibrium,  but  if  a lower  R.H.  were  specified  then  evaporation  would 

be  in  process  before  the  droplet  was  laser-irradiated. 

For  the  solution  droplets,  however,  this  drawoack  does  not  exist. 

- 13 

Considering  a 10  gm  particle  of  NaCl  in  a saturated  solution  (~  26% 

NaCl  by  mass,  molality  ~ 6),  one  obtains  p 1 = 1.  194  gm/cc  (also  a 

weak  function  of  concentration),  R = 0.425  pirn,  and  thus  from  Eq.  (59) 

P'r^Poo  ‘*9  °*  80’  so  that  such  a droplet  would  be  in  equilibrium  at  RH  = 80%. 

At  higher  humidity  the  droplet  would  increase  in  radius  and  decrease  in 

salt  concentration  until  it  established  a new  equilibrium.  Such  curves  are 

shown  in  Fig.  34,  the  shape  of  the  radius  vs.  RH  being  very  close  to  the 

4 

empirical  expression  given  by  Hodges  . This  is  shown  explicitly  in  Fig.  35 

where  the  predicted  droplet  radius  vs.  relative  humidity  as  determined  by 

the  variation  of  activity  with  molality  is  displayed  for  a NaCl  particle 
- 1 3 

of  10  gm.  The  droplet  size  shown  is  normalized  by  the  droplet  size  at 
saturation  (which  occurs  at  R.H.  **  80%).  Larger  salt  particles  would 
exhibit  the  same  growth  law,  however,  as  the  salt  particle  size  is  decreased 
surface  tension  effects  can  become  important  and  these  would  modify  the 

A 

growth  law.  Shown  for  comparison  is  the  growth  law  of  Hodges  , Eq.  (2), 
normalized  by  the  value  of  F at  80%  R.H.  The  comparison  is  quite  good. 
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Equilibrium  Radii  of  Sal,  Solution  Droplets  as  a Function  of 
Relative  Humidity  or  Supersaturation  (note  scale  change  above 
■ I00/o)-  Curws  are  reproduced  from  Ref.  10. 


NaCI  SOLUTION  DROPLET 
SATURATION  AT  R.H  =79.5% 


I-  R.  H./IOO 


Growth  Curve  for  Equilibrium  Radius  of  Salt  Solution  Droplet 
as  a Function  of  Relative  Humidity.  Solid  curve  is  calculated 
as  in  Fig.  34,  valid  for  NaCI  mass  > 10"^  g;  dashed  curve  is 
taken  from  Ref„  4. 


Unfortunately,  the  present  analysis  cannot  be  used  at  relative  humidities 
below  the  saturation  limit.  The  determination  of  p'  /p  is  equally  direct 
or  a mixture  of  solutes,  as  Chen  has  pointed  out  that  the  water  activity 
can  still  be  established  from  the  individual  solvent- solute  values  of  a 
with  appropriate  weighting  over  the  molality  of  each  solute  component! 

The  major  effect  of  inclusion  of  solution  properties  into  the 
modeling  of  aerosol  heating  would  be  to  reduce  the  vapor  pressure  at  the 
droplet  surface,  making  the  diffusion  boundary  condition  a function  of  both 
surface  temperature  and  salt  concentration  in  the  droplet.  As  vaporization 
proceeded,  the  solute  concentration  would  increase,  and  the  vapor  pressure 
would  become  an  even  smaller  fraction  of  the  saturated  vapor  pressure. 

-S°lute  Concentration  Gradients  within  the  Droplet 
The  evaporation  of  the  solvent  (H^O)  occurs  at  the  droplet  surface, 
so  that  solute  must  diffuse  towards  the  center  of  the  droplet  to  restore 
uniformity  of  concentration.  In  practice  this  will  almost  certainly  not 
occur  on  a sufficiently  fast  time  scale,  and  it  will  be  necessary  to  compute 
or  approximate  the  concentration  gradients.  A rough  estimate  of  the 
behavior  can  be  obtained  as  follows:  consider  a droplet  of  radius  R , with 
solute  concentration  Cq  gm/cc  at  the  surface.  If  the  solution  density  is 
PD  gm/cc  at  the  surface,  then  the  h^O  concentration  there  is  (p  - c ) 
gm/cc.  As  water  leaves  the  surface,  one  can  imagine  it  being  replaced 
by  water  flowing  from  the  center  of  the  droplet.  This  is  not  quite  the 
physical  situation  of  the  evaporating  droplet,  and  to  make  a steady  state 
situation  one  must  postulate  a source  flow  of  h^O  at  the  droplet  center. 

The  solute  then  has  a positive  gradient,  and  attempts  to  diffuse  back  to 
the  center  against  the  convective  velocity  of  the  water.  The  solute  con- 
centration at  radius  r is  then  given  by  In  (c/c  ) = A (1/R  - l/r),  and  the 


HO  flux  is  4 IT  R2  D n 
2 D 6 r 


^ ^ gm/sec.  The  constant  A 


may  be  established  by  specifying  the  radius  r'  at  which  c/c  is  some 
convenient  fraction,  say,  an  e-folding  distance,  so  that  -1  = A (1/R  - 
l/r'),  and  the  flux  = 4ITDPDR/(R/r'-l)gm/sec.  Now,  as  can  be  seen 
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from  Eq.  (34),  one  can  equate  the  H^O  mass  loss  to  0 4H  R3i  a /3A  H g/sec, 
where  0 is  the  fraction  of  absorbed  laser  energy  which  goes  into  vapori- 
zation. Equating  the  two  fluxes  results  in 


n Dp 

— - 1 = v D 

r 1 7 ’ 

PR  ia 
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Inserting  D = 0.  6 x l<f  cm  /sec  and  p D = U g/cc  for  NaCl  solutions. 


. 0123 
0 R2  I a 


For  small  values  of  laser  flux  R/r'  » 1,  which  implies  that  concentration 
gradients  are  very  small.  But  for  R2  I a * 1 We  have  R/r'  = 1 + ' 0123 
The  gradients  within  the  droplet  will  thus  depend  critically  on  the  fraction 
of  energy  going  into  vaporization;  if  B ~ .1,  say,  then  the  e-folding 

location  will  be  very  close  to  the  droplet  surface.  Under  such  circum- 
stances a real  droplet  would  develop  very  high  solute  concentrations  near 
the  surface,  probably  reaching  saturation  level.  Even  supersaturation  is 
a possibility,  since  the  recrystalization  process  would  require  nucleation 
and  could  well  be  kinetically  limited. 

3.  Ebullition 

The  possibility  of  nonuniform  heating  as  well  as  the  relatively  large 
temperature  differences  which  can  exist  within  the  droplet,  i.  e.  Eq.  (43), 
dictate  that  the  possibility  of  boiling  within  a droplet  be  considered.  If  a 
bubble  of  vapor  were  to  exist  within  a droplet,  instantaneously  in  mechanical 
equilibrium,  and  also  in  thermal  equilibrium  with  the  parent  liquid,  then 
the  pressure  within  the  bubble  would  be  p*  = P + 2a  /r*,  where  P is  the 
hydrostatic  droplet  pressure.  Taking  a = 58.  9 dynes/cm  at  T ~ 100°  C, 
results  in 
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P*/p  = t + 1.  178  (10"4/r«) 


(66) 


Thus,  even  a bubble  of  radius  1.2  Urn  would  require  a saturation  vapor 
pressure  of  2 atm  to  sustain  itself,  and  this  would  require  a liquid 

temperature  of  ~ 121  C.  A smaller  bubble  in  water  at  that  temperature 

would  collapse,  while  a larger  one  would  grow.  The  size  for  a self- 
sustaining  bubble  is  so  large  that  it  seems  extremely  unlikely  that  such 
bubbles  could  spontaneously  apnear  within  a 20-50  [lm  droplet  with  no 
solid  present  to  act  as  a nucleation  site.  The  same  statement  can  be  made 
even  more  strongly  for  a solution  droplet  with  no  insoluble  nucleus,  since 
the  solute  elevates  the  boiling  point,  and  raises  the  requirement  for 
superheat. 

However,  if  a supersaturated  solution  began  to  crystallize,  or 
if  an  insoluble  nucleus  were  present  in  the  droplet  it  is  clear  that  nonuniform 
heating  could  lead  to  internal  temperatures  which  would  yield  vapor  bubbles 
inside  the  droplet.  Indeed,  an  insoluble  nucleus  with  high  absorption 
coeiucient  could  enhance  the  nonuniform  heating  within  the  droplet  to  the 
extent  depicted  in  Fig.  32.  Thus,  this  aspect  of  the  ebullition  question 
merits  further  study. 

D.  Summary  of  Droplet  Effects 

All  the  effects  discussed  which  relate  to  surface  vaporization  of 
the  droplet  can  contribute  to  reducing  the  mass  loss  and  thereby  increasing 
the  transfer  of  heat  to  the  surrounding  air  by  conduction.  These  effects 
include  use  of  an  evaporation  coefficient  0.^  = . 033  instead  of  a » 1,  and 

lowering  of  vapor  pressure  at  the  surface  of  solution  droplets.  VFocusing 
of  the  laser  radiation  and  nonuniform  heating  effects  tend  to  increase  the 
heat  transfer  to  the  droplets  relative  to  the  volumetric  absorption  rates, 

and  to  lead  to  possibilities  of  ebullition  which  could  in  turn  lead  to  droplet 
shattering. 
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VI.  BEAM  PROPAGATION  EFFECTS 


This  section  is  concerned  with  the  resultant  effects  of  aerosol 
heating  and  vaporization  on  laser  beam  quality.  "Blooming"  of  the 
laser  beam  can  be  caused  by  two  separate  effects.  The  first  of  these 
is  that  the  heated  air  around  the  aerosol  will  have  a different  index  of 
refraction  from  ambient.  The  changes  in  the  index  of  refraction  across 
the  beam  width  can  be  a complicated  function  of  both  laser  intensity  and 
the  aerosol  size  distribution  function.  In  general,  the  phase  shift 
produced  in  a wave  propagating  through  this  nonuniform  medium  will 
depend  upon  these  variables  and  thus  there  will  be  variations  in  beam 
refraction  across  the  laser  beam  width.  Along  with  the  temperature 
variation  around  the  aerosols  there  will  also  be  a variation  of  water 
vapor  concentration.  Although  the  amount  of  water  vapor  vaporized 
from  the  aerosols  will  be  insignificant  compared  to  the  ambient  vapor 
concentration  the  former  will  be  sharply  peaked  around  the  aerosols  and 
in  that  region  can  provide  a major  part  of  the  total  index  of  refraction. 
While  the  index  of  refraction  of  water  vapor  is  different  from  that  of  air 
there  do  not  appear  to  be  any  measurements  nR  Q for  infrared  wave- 
lengths. It  is  clear,  however,  that  n^Q  would  have  to  be  significantly 
different  from  n^  for  the  water  vipor  perturbations  to  be  as  important 
as  the  temperature  perturbations  (note  at  X = 0.  8 Mm  the  ratio  (nf]  - 1)/ 

(nair  ” ~ °-  9)’  °nly  the  effect  of  temperature/density  variation^  will 

be  considered  in  the  following  analysis. 

The  prediction  of  beam  propagation  through  the  temperature  fields 
caused  by  heated,  vaporizing  aerosols  is  a very  complex  problem  which 
will  not  be  treated  in  detail  in  this  section.  However,  crude  estimates  of 
the  resulting  phase  shifts  will  be  made.  Both  Cook  and  Butts2  and  Lencioni 
and  Kleiman  have  examined  beam  propagation  through  the  temperature 
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fields  resulting  from  the  heating  of  nonvaporizing  aerosols.  Cook  and 
Butts  applied  the  first  Born  approximation  to  the  wave  equation  to  predict 
the  average  scattering  cross  section  of  an  aerosol  particle.  They  concluded 
that  scattering  from  the  hot  bubbles  would  only  be  important  at  large 
irradiation  times  (t  > 1 sec.  for  I = 10^  W/cm^). 

Lencioni  and  Kleiman  performed  a similar  analysis  for  the 
scattering  cross  section  but  also  predicted  the  average  phase  shift  a wave 
experiences  in  propagating  through  the  hot  air  bubble  about  a heated 
nonvaporizing  aerosol.  In  this  analysis  they  used  the  analytical  heat 
conduction  solution  and  included  beam  shape  effects.  They  made  the 
interesting  point  that  there  is  a hydrodynamic  sphere  outside  the  heated 
sphere  which  has  a positive  density  gradient  resulting  from  the  gas  which 
had  been  convected  outside  of  the  heated  bubble  at  the  speed  of  sound.  The 
inner  radius  of  this  hydrodynamic  sphere  occurs  somewhere  beyond 

r n/~4  kAt/p  Cp^,  the  point  where  the  conduction  solution  falls  off,  and 
the  outer  radius  is,  of  course,  given  by  ct. 

Their  main  conclusion  is  that  to  first  order  for  short  times  the 
average  phase  change  a wave  undergoes  in  traversing  the  hydrodynamic 
bubble  is  equal  and  opposite  to  that  experienced  within  the  heated  bubble. 

This  results  because  of  mass  conservation,  i.  e.  the  increase  in  mass  above 

ambient  in  the  hydrodynamic  bubble  must  balance  the  decrease  in  mass  in 
the  heated  bubble. 

A complication  in  this  is  that  the  hydrodynamic  bubble  is  not  well 
defined.  For  the  continental  distribution  described  in  Section  II  the  average 
spacing  between  aerosols  is  ~ 2 mm.  The  heated  bubbles  at  times  of  msecs 
are  typically  a few  hundred  Pm  in  radius  and  thus  will  not  overlap;  however, 
at  a time  of  1 msec  ihe  outer  radius  of  the  hydrodynamic  bubbles  will  be 
30  cm  which  is  of  the  order  of  the  laser  beam  width.  Thus,  not  only  will 
there  be  a strong  overlap  between  the  hydrodynamic  bubbles  of  neighboring 

aerosols  but  also  much  of  the  excess  mass  will  have  passed  outside  of  the 
beam  extent. 
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Because  of  this  one  could  conclude  that  the  phase  shift  encountered 
on  passing  through  the  heated  bubble  should  provide  a reasonable,  albeit 
somewhat  crude  estimate  of  the  total  phase  shift  encountered  by  the 
wave  after  an  irradiation  time  > , msec3.  In  any  event  in  the  case  of 

a vaporizing  molecule  mass  is  added  to  the  gas  at  the  expense  of  the  droplet 
and  the  constraint  of  gaseous  mass  conservation  does  not  apply  even  at 
short  times. 

Predictions  of  the  phase  shift  induced  in  a wave  propagating 
through  a heated  air  bubble  can  be  readily  performed.  The  index  of 
refraction  near  the  droplet  is  changed  due  to  the  change  in  air  temperature, 
ih.s  change  in  index  of  refraction  will  cause  a position  dependent  phase 
shift  for  rays  passing  through  the  heated  air.  At  constant  pressure 
the  index  of  refraction,  n,  at  temperature  T is  related  to  the  ambient 
index  of  refraction,  n,x  , at  temperature  Ta  by: 


n - i 

n - 1 

00 


(67) 


The  present  analysis  will  be  restricted  to  cases  where  T - Tro  is  much  less 

than  T„.  Then  expanding  Eq.  (67)  in  a binomial  series  and  keeping  only 
the  first  two  terms  results  in 


(n^  - 1 ) A T 
n = n 

00  T (68) 

00 

where  A T = T - . This  index  of  refraction  will  be  used  to  calculate 

tne  optical  path  through  a heated  air  bubble. 

I he  optical  path  length  of  a ray  is  defined  by 

OP  = f nds  (69) 

C 
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where  the  ray  is  taken  to  follow  a curve  C through  the  medium  and  ds 
is  the  differential  element  along  this  curve  C.  Given  an  expression  for  the 
temperature  field  about  the  droplet  Eq.  (68)  and  (69)  can  be  combined  to 
evaluate  OP.  As  discussed  in  Section  IV  the  radial  temperature  distribu- 
tion about  a heated  aerosol  droplet  at  momra,e  laser  intensities  such  that 
convection  cculd  be  neglected  was  given  by 


LT  (r) 


£ T (R)  R 
r 


(25) 


foi  the  legime  R < r < rQ,  where  R is  the  aerosol  droplet  radius,  r is 
the  radial  distance  from  the  droplet  center,  and  rQ  is  a characteristic  heat 

p Combining  Eq.  (25)  with  (68)  results 

Pi 


conduction  length,  V 4 k.t/o  C 


in 


<n„  - U Or  (R)  R 


Tm  r 
00 


R < r < 


r > 


(70) 


= 1.  36 


r •:  R 


In  general  in  calculating  the  total  path  length  the  shape  of  the  curve 
must  also  be  evaluated.  In  the  following  analysis  C will  be  taken  as  a strai 
line  and  it  will  be  seen  that  the  increase  in  path  length  over  ambient,  A OP, 
is  sufficiently  small  to  justify  this. 

A heated  air  bubble  of  radius  r^  is  sketched  below.  If  a ray 


C 

ght 


- 


enters  at  point  and  exits  at  point  it  can  readily  be  shown  from 
Eqs.  (69)  and  (70)  that 


2 (n  - 1)  AT  (R)  R 


OP  n (z  - z.)  + -- 
00  2 1 


In  [ ta n ( 6/2  ) ] (71) 


providing  that  x;  > R.  (Entry  points  P^  corresponding  to  x^  < R 
will  be  excluded  from  the  analysis).  Thus 


2 (n  - 1)  A T (R)  R 


A OP  = 


In  [ tan  (9  /2)  J . 


(72) 


Note  that  the  quantity  A OP  is  not  independent  of  bubble  size  since 

9 = sin  1 (x/r  ). 

1 o 

A detailed  phase  shift  prediction  would  require  tracing  a wave 

propagating  through  an  ensemble  of  such  heated  bubbles,  as  expected 

from  a typical  aerosol  distribution,  including  the  wave  refraction  and 

random  entry  sites,  , at  each  bubble.  Instead  the  average  phase 

shift  encountered  in  propagating  through  a bubble  will  be  calculated 

as  an  estimate  of  the  actual  solution. 

The  quantity  In  [tan  (0  / Z)~]  can  vary  considerably  for  values 

of  x /r  between  R/r  and  unity.  However,  it  can  be  shown  that  the 
1 o o 

average  value  of  A OP  defined  by 


A CP 


r 

o 

A OP  2 TT  x dx, 
R 1 1 


2 2 
TT  (r  - R ) 
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(73) 
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is  given  by 


A OP  *» 


2 (n  - 1)  AT  (R)  R 


(74) 


So  far  this  analysis  is  similar  to  that  of  Lencioni  and  Kleiman^;  however, 
in  the  case  of  a water  dr  plet  the  quantity  A T (R)  will  no  longer  be 
linearly  proportional  to  la  for  values  of  R la  > 10  - 1 0 cal/cm  - sec. 

As  a specific  example  of  the  average  phase  shift  encountered  by 
a wave  propagating  through  an  atmospheric  aerosol  distribution  consider 
the  case  of  the  standard  continental  distribution,  Eq.  (1),  at  an  R.H.  of 
q8?o.  The  average  phase  shift  is  given  by 


R, 


<0  = Z TT 


2 2 TT 


2 TT  f 

1 4 


A OP 


d N 
A 

d R 


dR 


(75) 


where  z is  the  range  in  cm  and  dN  /dR  is  the  aerosol  distribution  given 
by 


d N _ i l 

A . K 3,5x  10  .4 

d R - 98^u)  5 > part/cm  (76) 

R 

for  0.5  Um  < R < 100  [im. 

From  reference  to  Fig.  4 it  can  be  seen  that 


A T (R) 


R2  la  „ 2 _ _ _ -2 

^ ^ > R Id  < 4.5x10  cal/cm-sec 


80  K;  R la  > 4.  5 x 10  cal/cm-sec 


(77) 
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or  approximately  1/3  of  the  heat  conduction  solution  for  R I a < 1.5  x IQ'2 

cal-cm-sec.  Considering  as  a specific  example  I = 104  watts/cm2 

and  a = 113cm  , then  R la  = 4.5x10  cal/cm-sec  at  R = 4 Mm.  Therefore, 

coupling  Eqs.  (74)  - (76)  and  the  definition  that  r 2 = 4 k t/pCn  results  in 

o A P A 


<o  = 


1 6 tt  z k t (n  - 1 
A ® 


PCp  X T» 
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4x10 


-4 
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3.  4 x 10  1 *Ia 


0.  5x10 
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9 k R 
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3.  5x  1 0"  1 Lx  80 


dR 
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(78) 


For  this  particular  case  the  dominant  contribution  is  provided  by  the 
first  integral,  for  which  0.  5 < R < 4 m,  and  the  result  may  be 

approximated  as 


<o 


-9 

1.3x10  zt(n  - 1)  la 


P C, 


X T 


(79) 


53  . 

Inserting  la  = 2.  7 x 10  cal/cm  -sec,  no  - 1 = 3 x 10  and  X = 3.  8 |im 
yields 


(<ty/z  w 300  t km  (g( 

corresponding  to  phase  shifts  of  order  unity  at  times  of  msecs  and 
ranges  of  ~ 1 km. 

As  long  as  the  only  droplets  which  can  reach  the  boiling  point  come 
from  the  tail  of  the  aerosol  distribution  function,  the  above  result  (Eq.  (79)  ) 
is  not  strongly  affected  by  the  uncertaii.ties  discussed  in  previous  sections 
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Of  this  report.  For  instance,  use  of  a vaporization  coefficient  of 

3 v = ' 033  WOuld  raise  the  of  A T (R)  in  Eq.  (77)  by  at  most  a 

factor  of  3.  On  the  other  hand,  at  higher  laser  intensities  the 
vaporization  limit  is  approached  at  smaller  droplet  radii,  and  the 
details  of  the  boiling  phenomenon  become  more  important  in  the 
determination  of  beam  quality. 

It  should  be  emphasized  that  there  is  no  justification  for  adding 

and  averaging  the  phase  shifts  in  the  above  analysis.  This  was  done 

only  in  order  to  provide  a gross  estimate  of  the  effect  of  aerosol  heatir.j 

on  beam  quality.  More  realistic  predictions  must  await  a more  detailec 
analysis. 

Furthermore,  phase  shifts  resulting  from  molecular  absorption 
can  be  larger,  particularly  if  the  H^O  vapor  continuum  absorption  is  as 
strong  as  it  appears  to  be.  However,  as  pointed  out  in  Section  II  aeroso 
mass  loadings  higher  than  the  standard  continental  distribution  are 
commonly  observed  and  can  be  the  dominant  absorbers.  It  should  be 
emphasized  that  the  functional  form  of  Eq.  ( 79)  is  similar  to  that 
resulting  from  uniform  absorption  because  the  individual  phase  shifts 

were  averaged.  This  is  functionally  the  same  as  averaging  the  density 
change  in  the  gas. 
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VII.  SUMMARY 


A detailed  model  has  been  developed  to  predicL  the  temperature 
and  vapor  fields  resulting  from  the  laser  irradiation  of  small  water 
droplets.  The  water  droplets  are  considered  to  be  representative  of 
hygroscopic  aerosols  at  high  relative  humidities.  The  model  includes 
the  effects  of  vaporization,  heat  conduction  and  convection.  [’he  model 
was  developed  under  the  assumption  of  a constant  pressure  gas  and  thus 
is  valid  only  up  to  tne  isobaric  vaporization  limit. 

Both  steady  state  and  time  dependent  solutions  for  the  model  have 
been  discussed  in  some  detail.  It  was  shown  that  the  predicted  profiles 
of  increase  in  gas  temperature  and  water  vapor  scale  inversely  with 
distance  from  the  droplet  for  moderate  laser  intensities  and  particle  sizes. 

F’or  larger  laser  intensities  and/or  particle  sizes  convection  began  to 
dominate  conduction  and  diffusion  resulting  in  more  bowed  profiles.  It 
was  demonstrated  that  the  model  was  specifically  valid  only  if  the  evaporation 
coefficient  for  water  a approached  unity.  It  was  pointed  out  that  recent 
evidence  suggests  tliata  ^ is  considerably  smaller  than  unity  and  if  this  is 
so,  the  model  must  he  modified  accordingly. 

A number  of  effects  corresponding  to  real  aerosol  properties  not 
included  in  the  model  were  discussed  in  some  detail.  These  include  nonuniform 
absorption  and  heating,  solution  effects  including  recrystallization  and 
ebullition.  It  is  concluded  such  phenomena  could  be  important  under  extren  e 
conditions  corresponding  to  heavy  vaporization  and/or  droplet  superheating. 
LsLin; ates  of  the  characteristic  phase  shifts  a wave  will  encounter  in 
propagating  through  a heated  aerosol  distribution  were  also  presented. 

There  are  several  areas  of  this  problem  which  warrant  further 

study: 

a)  The  estimates  of  characteristic  phase  shifts  are  sufficiently 

large  to  require  that  a more  accurate  prediction  be  performed. 
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b)  Model  predictions  should  be  made  for  the  case  of  a 0.33. 
Furthermore,  a firmly  established  value  of  this  quantity 
would  be  valuable. 

c)  The  absorption  coefficient  of  liquid  water  near  3.  S pm 
should  be  re-examined  along  with  its  temperature  dependence 
in  order  to  clear  up  present  uncertainties. 

d)  Lastly,  the  water  vapor  continuum  absorption  coefficient 
and  index  of  refraction  near  3.8  pm  should  be  measured 
near  room  temperature. 
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